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ill  ^ 
PREFACE 


This  Report  Is  one  in  a  series  entitled  "Effects  of  Strong  Eicplosions" 
and  intended  to  provide  an  overview  of  current  Soviet  research  activities 
in  that  field.  It  falls  within  the  scope  of  a  continuing  program,  spon¬ 
sored  by  the  Defense  Advanced  Research  Projects  Agency,  which  undertakes 
the  systematic  coverage  of  selected  areas  of  Soviet  scientific  and  tech¬ 
nological  literature. 

The  effects  of  strong  explosions,  particularly  nuclear  explosions, 
and  their  simulation  under  laboratory  and  fi'Sld  conditions  involve  a  very 
broad  range  of  phenomena  and  parameters  such  as  high  temperatures,  pres¬ 
sures,  electromagnetic  energy  densities,  interaction  of  energy  with 
materlale,  shock  waves  in  gases,  liquid,  and  solids,  and,  in  general, 
events  occurring  under  conditions  of  a  high  degree  of  mechanical,  thermal, 
and  radiation  stress.  Tlie  series  of  Rand  Reports  of  which  this  is  the 
fourth  to  be  published  will  cover  a  variety  of  the  pertinent  topics  as 
they  are  reflected  in  the  currently  available  Soviet  literature,  by 
supplying  abstracts  from  and  summaries  of  the  scientific  and  technolog¬ 
ical  publications  of  the  USSR. 

The  material  in  the  present  Report  is  derived  exclusively  from  Soviet 
technical  writings  that,  with  a  few  exceptions,  appeared  in  1971  and  the 
early  part  of  1972.  The  abstracts  here  Included  were  prepared  by  Informa¬ 
tics,  Inc.  They  have  been  arranged  by  the  compiler  according  to  fairly 
general  subject  areas,  and  vary  in  length  depending  on  the  pertinence 
and  significance  of  the  individual  papers. 

The  information  assembled  in  this  Report  covers  our  most  recent 
knowledge  of  Soviet  research  on  the  effects  of  explosions,  and  should 
prove  useful  to  American  scientists  working  in  this  and  related  fields 
of  research. 
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I.  NUCLEAR  EXPLOSION  EFFECTS 


Kukhtevich,  V.  I. ,  I,  V.  Goryachev,  and  L.  A., 

Trykov,  Zaebchita  ot  pronikayushchey  radiatsii 
yadernogo  vzryva  (Protection  against  the  penetrating 
radiation  of  a  nuclear  explosion  ) .  Moscow,  Atomizdat, 

1970.  190  p. 

An  engineering  design  method  is  presented  for  the  protection 
of  ground  structures  against  the  penetrating  radiation  of  a  nuclear  explosion. 

The  book  is  based  on  Soviet  and  Western  open  source  materials  published  in 
1958-1969;  Soviet  sources  make  up  about  50%  of  the  123  references. 

Physical  aspects  of  nuclear  explosion  are  treated  only  briefly, 
since  they  are  covered  well  enough  in  several  recent  Soviet  and  Western  mono¬ 
graphs.  Emphasis  is  placed  on  neutrons  and  gamma-ruys  produced  at  the 
atmosphere  -  ground  boundary  and  their  penetration  through  a  typical  protective 
shield.  Numerous  published  data  on  characteristics  of  ionizing  radiation  are 
assembled  for  use  in  the  design  of  protection  against  nuclear  radiation. 

Formulas  and  numerical  data  are  given  for  calculation  of  the 
dose  rate  cf  primary  gamma  radiation  propagating  from  a  point  source  and 
from  a  spherical  radioactive  cloud  through  an  infinite  atmospheric  medium. 

The  effect  of  a  shock  wave  created  by  an  explosion  on  the  propagation  of 
primary  gamma-rays  is  evaluated  in  the  cases  of  aerial  and  ground  explosions. 
Formulas  are  also  given  for  the  dose  rates  of  primary  and  secondary  gamma 
-radiation  with  allowance  for  the  shock  wave  effect,  as  well  as  the  time, 
explosive  force,  distance  from  the  explosion  center,  and  ascent  velocity  of 
the  radioactive  cloud.  Numerical  values  for  an  integral  dose  of  gamma 
radiation  at  a  given  distance  from  the  explosion  center  are  described  by  formulas 
for  low  and  high  explosive  forces,  and  plots  of  gamma-ray  accumulation  rate 
versus  time  are  shown.  Theoretical  data  are  included  on:  spatial  distribution 
of  slow  neutrons  released  by  a  nuclear  explosion  and  formed  in  air  by  moder¬ 
ation  and  diffusion  of  deUyed  neutrons;  integral  neutron  albedo  from  the 
ground  surface  and  coefficients  of  the  angular  distribution  function  of  scattered 
neutrons;  and  characteristics  of  neutrons  reflected  from  the  ground  by  a 
nuclear  explosion  in  the  atmosphere.  Calculated  data  are  also  presented 
on  integral  gamma -radiation  albedo  in  a  nuclear  explosion. 

Other  tables  give  theoretical  and  experimental  data  on  transmission 
of  nonscattered  neutrons  and  total  neutron  dose  through  an  iron  shield  of 
varying  thicknesses  and  at  different  angles  of  incidence;  transmission  of  a 
total  dose  of  neutrons  through  water,  polyethylene,  and  iron  shields;  trans¬ 
mission  of  nonscattered  neutrons  from  a  nuclear  reactor  through  an  iron, 
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concrete,  or  polyethylene  shield;  and  angular  distribution  of  the  dose  rate 
of  scattered  neutrons  beyond  an  iron  or  polyethylene  shield.  Experimental 
and  theoretical  data  are  plotted  and  tabulated  on  angular  distribution  (differ¬ 
ential  albedo)  and  integral  albedo  of  neutrons  back- scattered  from  water, 
ground,  and  iron  at  different  neutron  beam  angles  of  incidence. 

Theoretical  mass  coefficients  of  attenuation  of  nonscattered 
gamma- radiation  by  different  materials  are  tabulated  and  discussed.  Theo¬ 
retical  transmission  coefficients  of  gamma- radiation  scattered  by  iron  shields 
of  varying  thicknesses  and  at  different  angles  and  characteristic  angles  of 
scattering  for  specific  energies  and  materials  are  tabulated.  Formulas  used 
to  calculate  the  tabulated  data  are  given.  Empirical  formulas  .ire  cited  for 
calculating  integral  energy  albedo  of  gamma  quanta  back- scattered  from 
protective  shields.  The  approximate  formula  for  bremsstrahlung  energy  of 
beta-electrons  is  given  and  data  on  the  relative  contribution  of  spectral 
components  of  bremsstrahlung  were  tabulated.  Secondary  gamma- radiation 
yield  beyond  the  shield  is  analyzed,  and  the  numerical  values  of  the  constants 
for  calculation  of  the  yield  are  tabulated  for  water,  iron,  and  concrete. 
Calculated  yield  functions  for  different  shield  thicknesses  are  plotted. 

Two  methods  of  designing  the  protection  of  ground  structures 
are  outlined:  (1)  an  exact  method  involving  the  use  of  computers  and,  (2)  a 
simplified  method  using  approximate  formulas.  Both  methods  require  a 
knowledge  of  the  barrier  and  geometry  factors  of  radiation  attenuation.  The 
formulas  for  calculation  of  these  factors  are  given,  and  the  sequence  of 
mathematical  operations  associated  with  the  two  methods  is  described.  The 
precise  method  consists  of  six  consecutive  steps. 

The  effects  are  evaluated  of  the  multiple  scattering  of  radiation 
on  the  inner  walls  of  a  structure. 

Appendices  contain  8  tables  of  barrier  factors  of  transmission 
of  neutrons  and  gamma-radiation  through  various  materials,  and  two  tables 
for  yields  of  trapped  gamma -radiation  for  various  materials. 


Kozlov,  V.  I.,  and  Yu.  G,  Shafer.  On  a  possible 
mechanism  for  the  localization  of  fission  products 
from  a  high-altitude  nuclear  explosion.  Kosmicheskiye 
issledovaniya,  no.  4,  1971,  630-631. 

Recorded  particle  data  on  the  Cosmos  3,  5  and  6  satellites 
following  the  U.  S.  Starfish  high-altitude  detonation  in  July,  1962  showed  a 
concentration  of  fission  products  at  the  conjugate  point  to  Johnson  Island. 
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This  concentration  was  detected  at  high  altitudes  and  persisted  for  an  appreciable 
length  of  time.  In  the  present  paper  the  authors  use  the  Starfish  results  for  a 
model  with  which  they  analyze  the  type  of  trapped  particles,  their  transport 
oharacteristics,  and  the  geometry  of  the  trap  region  at  the  conjugate  point, 

A  magnetic- gravitational  trap  is  postulated  in  the  conjugate  point 
region  which  will  act  as  a  potential  well  for  the  heavier  charged  particles.  At 
some  radial  distance  rQ  from  earth  center  the  combination  of  gravity  and  repulsion 
potentials  will  have  a  minimxim,  which  determines  the  depth  of  the  potential 
well;  in  the  proposed  model  this  is  shown  to  be  E  j_  =  mg  r/3,  where  m  -  ion 
mass.  Assuming  particles  with  atomic  weight  of  100  and  an  altitude  of  400  km, 
we  find  that  E  ,  ^  30  ev,  i.  e.  ,  the  portion  of  the  dispersing  cloud  of  fission 

products  having  energies  on  the  order  of  30  ev  will  be  trapped  a.t  the  400  km  le\el. 
This  level  is  consistent  with  the  values  obtained  by  Neff  et  al  [ij  and  Colgate  [2] 
at  the  conjugate  point. 

■  The  arrival  time  of  the  major  portion  of  charged  particles  is  generally 
a  function  of  plasma  dispersion  in  the  geomagnetic  disruption  at  the  detonation 
point;  Askar'yanf33  gives  the  velocity  as  on  the  oi'dei*  of  cm/sec  for  the 
Starfish  test,  which  for  a  force  line  length  of  about  4  x  10  cm  yields  an  arrival 
time  of  40  seconds. 

The  authors  next  examine  the  factors  determining  the  duration  of  the 
trap  region.  As  the  analysis  shows,  the  size  of  the  well  is  a  direct  function  of 
particle  charge,  hence  the  well  ass\imes  a  "funnel"  configuration  according  to 
trapped  charge  level.  The  confinement  time  is  given  approximately  by 


’■j,  =  0.785  X  In  l^max  , 

B 

where  is  mean  ion  relaxation  time  and  B  is  magnetic  inductance,  on  tlie 
assumption  that  the  main  loss  mechanism  is  coulomb  scattcring  in  the  funnel  region. 
If  a  Maxwellian  temperature  distribution  of  atmospheric  ions  is  assumed,  then 
a  numerical  solution  for  under  the  assumed  atomic  weight  and  ionization  level 
falls  in  the  range  of  lO^.-lO^  sec  in  the  250--1000  km  altitude  range.  This  suggests 
that  the  cone  region  is  sustained  by  continuous  replenishment  of  particles  from 
the  upper  energy  regions. 

Finally,  it  is  concluded  that  dissipation  of  the  funnel  region  owing 
to  drift  of  trapped  particles  is  negligible  (0.2  -  0.  5  deg/day),  which  was 
corroborated  by  experimental  evidence  that  the  trap  region  acted  as  a  compact 

/3  and  7  radiation  source  for  a  substantial  period. 
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This  paper  is  the  first  extensive  analysis  to  appear  of  Soviet 

968.  At  the  time  of  the  explosion  two  Cosmos  satellites  were  in  orbit  and 
carrying  tyP®8  radiation  monitoring  equipment:  Cosmos-259  (apogee 

Cosmo8-262  (apogee  818,  .perigee  263  km).  The 
orbital  planes  of  both  were  inclirw^  at  49°  to  the  equatorial  plane.  The 

STrl’;  tv™  triple-coincidencn  panoramic  telescope, 

featuring  type  STS- 5  gas  discharge  counters  (identified  as  the  STS  channel),  and 
n  unshielded  STS- 5  counter,  identified  as  the  SSVS  channel,  mounted  on  the 
end  of  a  2  meter  boom.  The  panoramic  system  recorded  only  charged  particles 

fn*^«  d"  hT  ^  electrons  >  12.  5  MeV;  the  single  SSVS  counter  waL 

intended  for  protons  >  10  MeV  and  electrons  >  50  KeV;  it  could  also  register 

conluH  «  f  f  intensity  array 

br3m  All'^^  3°“*’  PK  ‘^i^^harge  counters  connected  in  parallel  and  shielded 

oy  3m  Al  and  3  mm  Pb  screens. 


i-.-  ..  .  Th„  equipmept  was  uaed  to  transmit  real-time  data  on  radiation 

thl  time’' th  d  USSR,  Mongoiia  and  China,  around 

!h!  Jn^  ^«o"»t>on.  For  the  purpose  of  illustration  the  authors  have  chosen 

tadiau™"  ’i’hi  d  r”  “•**’*'  neur-eart 

as  shown'in  Fin.  1  3  ’  P.”"'".*"*  according  to  satellite  trajectory 

.  coin'd  increase  in  count  rate  N  over  that  of  a  "quiet"  day 

i.e. ,  December  18  and  19,  1968.  ^ 

recorded  d.7-^*  4howthat  the  shielded  single  counter  (SSC)  in  the  Co8mos-259 
recorded  during  two  minutes  an -80%  increased  intensity  in  its  196th  orbit,  i.e.  , 
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about  lb  hours  after  the  explosion;  this  recording  was  made  at  an  altitude  of 
about  700  km  over  the  Mongolian  Peoples  Bepublic  in  an  area  whose  coordinates 
were  43-45°N.  Lat.  (geographic)  and  lOO-llO^E.  Long.  The  counting  rate  of 
tie  triple-coincidence  telescope  increased  at  the  same  timel  Twenty-four  hours 
after  the  explosion,  over  approximately  the  same  territory  in  the  309th  and  210th 
orbits  above  the  eastern  territory  of  China,  the  satellites  noted  an  increase  of 
intensity  (AN/N)ggQ  =  20  and  60%,  (AN/N)g.j.g  =  50  and  40%,  respectively.  In 
the  211th  orbit,  the  counting  rate  increments  of  the  single  shielded  counter 
reached  1.5-2  pulses  per  second,  which  amounts  to  30-40%  over  background. 

On  December  27  in  the  196th  orbit,  beginning  with  2113-2117  hours, 
the  external  SSVS  counter  rates  increased  by  a  factor  of  2-8;  a  maximal  increase 
by  a  factor  of  10  was  noted  at  2115  hours  30  seconds.  Subsequently  the  count 
rate  began  to  fall  off  and  at  2118  hours  only  background  oscillation  v.'as  recorded. 

In  the  197th  orbit,  the  count  rate  increased  by  a  factor  of  2-2.5.  On  December 
28  in  the  209th  orbit,  the  satellite  recorded  only  the  cosmic  ray  background  while 
in  the  210th  orbit,  whose  trajectory  is  near  the  196th  orbit,  in  a  real-time 
transmission  regime,  the  count  rate  increased  by  a  factor  of  1.5-2.  The  maximum 
"excess"  of  the  count  rate  in  the  211th  orbit,  whose  trajectory  almost  coincided 
with  the  197th  orbit,  reached  100%. 

Analogous  results  are  reported  for  the  Cosmos-262  count  rate,  e.  g.  , 
262 's  counters  registered  a  doubled  intensity  23  hours  after  the  e.xplosion;  the 
recording  was  also  made  over  the  Mongolian  Peoples  Republic  in  an  area  whose 
coordinates  were  48-47°N,  Lat.  and  95-1059E.  Long.  ,  at  an  altitude  of  270  km. 

It  is  noteworthy  that  no  increase  in  radiation  intensity  was  recorded 
when  the  Co8mos-259  satellite  passed  in  its  209th  orbit  over  the  magnetic 
conjugate  of  the  Chinese  explosion  point,  which  is  to  be  expected  in  low  altitude 
explosions.  Consequently,  there  was  no  direct  ejection  of  explosion  products 
into  the  force  tube.  A  further  theoretical  analysis  is  given  to  define  the  absorption 
mechanism  taking  place  in  tjje  vicinity  of  the  detonation  point;  the  fact  that 
neither  Cosmos  satellite  intercepted  fission  fragments  is  taken  as  evidence  of 
strong  atmospheric  absorption,  and  on  this  basis  it  is  concluded  that  detonation 
occurred  at  not  over  20  km  altitude. 

The  level  of  protons  above  190  MeV  registered  by  the  STS  counters 
suggests  that  there  may  have  been  some  release  of  them  from  the  radiation 
belts  owing  to  the  disturbing  effect  of  the  Chinese  detonation.  Using  calculations 
from  the  U.  S.  Starfish  test,  the  authors  verify  this  mechanism  as  a  possible 
additional  result  of  the  test. 


monuclear  explosion  according  to  data  from  the  panoramic 
triple-coincidence  telescopes  (STS)  on  board  Kosmos-259. 

I  —  196th  orbit,  12-27-1968;  II  —  210th  orbit,  12-28-1968; 
III  -  209th  orbit,  12-28-1968. 


counter  on  board  Kosmos-259.  I  —  197th  orbit,  12-27-1968; 
II  —  211th  orbit;  12-28-1968;  III--  196  th  orbit,  12-27-1968; 
IV  —  210th  orbit,  12-28-1968;  V  —  209th  orbit,  12-28-196a. 
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II.  EXPLOSIONS  IN  VARIOUS  MEDIA 


Rodionov,  V.  N. ,  and  A.  P.  Sukhotin. 
Parameters  of  elaatic  waves  generated  by 
spherical  explosion  in  metals  of  different 
strength.  FGiV,  no.  1,  1971,  142-146. 


An  experimental  study  was  made  of  the  effect  of  yield  point  a 
in  uniaxial  elongation  of  A-00  aluminum  (1),  D-16  annealed  (2),  and  D-16 
hardened  duralumin  (3)  on  the  parameters  of  an  elastic  compression  wave 
generated  by  a  0. 2  g  spherical  charge  of  TEN.  A  diagram  of  the  experiment 
is  shown  in  Fig.  1.  The  metal  samples  had  similar  mechanical  characteristics, 


Fig.  1.  Experimental  setup:  1  -  explosive 
charge;  2  -  metal  block;  3  -  pickup  plate; 

4  -  guard  ring;  5  -  insulators;  6  -  cable. 


with  the  exception  of  differences  in  strength.  The  measured  0  was  450, 

1600,  and  3«700  kg/cm  for  samples  (1),  (2)  and  (3).  Velocity  rise  time 
T  j,  positive  phase  duration  T2,  and  maximum  mass  velocity  U,„  were 
measured  on  oscilloscope  traces  of  free  metal  surface  displacement.  The 
experimental  plots  tj  (r),  t2  (r),  (Fig.  2),  and  lgU|^(lgr)  indicate  that 
elastic  waves  are  generated  by  the  explosion,  and  the  observed  differences 
In  temporal  wave  parameters  are  the  effect  of  strength  only.  This  effect 

was  also  evident  in  variations  between  the  measured  volume  of  the  explosion 
cavities  of  the  three  metal  samples  (1.3,  0.75,  and  0.42  cu  cm,  respectively). 
The  linearity  of  the  calculated  plots  of  the  relative  maximum  mass  velocity 
versus  distance  from  the  explosion  center  confirms  that  the  extent  of  plastic 

deformation  is  determined  by  the  static  a  of  the  metal.  The  motion 
pattern  and  state  of  the  metal  in  the  elastic  region  was  reconstructed 
from  the  experimental  data  and  the  equation  of  plasticity.  The  effect 
of  a  on  the  motion  of  the  boundary  line  R,|,  at  which  the  difference  of 
principal  stress  is  equal  to  a  of  the  metal  is  shown  in  Fig.  3.  An  increase 

in  a  affects  both  the  amplitude  and  temporal  parameters  of  R^i  motion. 


to  t,  M«cc 


Fig,  3.  The  effect  of  9  of  the  2  and  3 
metalD  on  boundary  line 
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Krupnikov,  K.  K. ,  V.  F.  Kuropatenko,  A.  T. 
Sapozhnikov,  B.  N.  Slmanov,  and  V.  A.  Simonenko. 
Calculation  of  explosions  in  a  medium  with  poly¬ 
morphic  phase  transition.  DAN  SSSR,  v.  202, 
no.  2,  1972,  300-301. 


An  explosion  in  a  quartz-like  medium  was  calculated  on  a 
computer  using  a  RAND  (rashet  adiabaticheskikh  nestatsionarnykh  dvizheniy) 
program  for  acnstationary  adiabatic  motion.  The  program  takes  into  account 
the  polymorphic  phase  transition  of  quartz  into  "stishovit".  A  two-liquid 
model  was  adopted  for  the  medium  with  a  gas  bubble  in  its  center;  high  pressure 
expansion  of  the  bubble  simulated  a  point  explosion.  In  the  region  of  continuous 
flow,  pressure  P  and  internal  energy  E  of  each  phase  were  calculated  from 
the  equations 

dF  +  PdV  — 0, 


(1) 

(2) 


where 


8  =  p/po  =  0.605. 


ni  =  4,  »i2  =  2,7,  Fi  =  r2  =  0.55,  Eq  =  const,  P,  V,  and  E  were  corre¬ 
lated  by  a  =  P(T)  equilibrium  curve  and  the  equation  of  state  of  the  mixture, 

jg-£.  =  (rdP(r)/dr-p(r)j(r-F,) 

where  V*  and  Eitr  arc,  respectively,  specific  volume  and  internal  energy 
of  one  of  the  phases  along  the  equilibrium  curve.  In  many  cases,  expansion 
shock  and  bifurcated  compression  waves  appear  at  the  phase  mixture  boundary. 
Allowance  was  made  in  the  RAND  program  for  these  wave.^.  The  calculated 
shock  wave  parameters  were  plotted  against  a  dimensionless  variable 
r  =  r/rg,  v  here  r^  =  [q/(  p  qCJ)]  (Fig,  1).  The  plots  show  that  P  and  U 


Fig.  1.  Dimensionless  pressure  P  and 
velocity  0  versus  time:  1  -  before  bifuration 
of  the  first  shock  wave;  2,  3,  and  4  -  after 
*  bifuration. 


behind  the  first  shock  wavefront  increase  significantly  before  dropping 
stepwise  at  the  front  of  the  expansion  shock  wave. 


Shurshalov,  L,  V.  Calculation  of  powerful  underwater 

explosions.  MZhiG,  no.  5,  19'’!,  36-40. 

A  problem  on  powerful  underwater  explosions  is  solved  by  finite 
-difference  equations,  with  the  introduction  of  an  artificial  viscosity  Q  =  Qj  +  Q2, 
where  Qj  and  Q2  are  linear  and  quadratic  components  of  Q.  In  contrast  to 
relatively  weak  explosions,  a  full  thermodynamic  description  of  the  properties 
of  water  over  a  broad  range  of  pressures,  temperatures  and  densities  is  required 
In  the  present  problem.  A  numerical  solution  is  developed  using  an  equation 
of  state  derived  earlier  by  the  author  (MZhlG,  no.  4,  1967,  184)  to  describe 
the  properties  of  water  at  pressures  over  1,000  atm.  At  a  pressure  below  1,000 
atm,  the  Talt  equation  of  state  was  used. 

A  spherical  envelope  (v  =  3)  containing  a  strongly  heated  compressed 
gas  and  located  in  an  infinite  space  filled  with  water,  with  specific  volume  V2 
at  a  pressure  pj,  served  as  a  mathematical  model  of  an  underwater  explosive 
charge.  Motion  of  water  and  gas  was  calculated  from  the  time  t  =  0,  when  the 
envelope  disappears  Instantaneously.  A  set  of  Lagranglan  equations  of  motion, 
which  included  the  term  Q,  and  ec^uations  of  state  of  detonation  products  and 
water  were  solved  by  the  finite>difference  method.  A  numerical  calculation 
using  the  finite-difference  equations  was  performed  for  explosion  of  a  spherical 
TNT  charge  of  =  1,  5  g/cm^  at  about  60  m  depth.  Characteristic  values 
selected  were:  p*  =  10^  kg/m2,  V*  =  10-2  x  sec2,  and  1  =  ro,  where 

to  is  the  initial  radius  of  the  charge.  The  coefficients  Cj  and  C2  in  the  expressions 
Ql  =  Cja  AU/V  and  Q2  =  AU^/V  (where  AU  =  U(R  +  AR)-UR,  U,  V,  and  R 
are  dimensionless  velocity,  specific  volume,  and  Euler  coordinate  respectively, 
and  a  is  the  local  sound  velocity)  were  assumed  to  be  Ci  =  0.  5  for  extremely 
powerful  or«0  for  weak  waves  and  C2  =  2.  The  calculated  hydrodynamic  para¬ 
meters  are  plotted  in  Fig.  1-3,  where  r=  t/t*,  dashed  lines  describe  variations 
of  the  shock  wave  parameters,  and  Rj  =  1  corresponds  to  the  boundary  of  a  gas 
bubble.  The  dots  and  circles  represent  the  data  calculated  with  a  time  step 
twice  (dots)  and  one  half  (circles)  that  used  in  calculation  of  data  represented  by 
the  curves.  Comparative  analysis  shows  that  the  calculations  are  natisfactorily 
accurate.  The  flow  parameters  were  also  calculated  for  dispersion  in  water  of 
a  spherical  (v  =  3),  cylindrical  (v  r  2),  or  a  plane-parallel  (v  =  1)  volume  of  a 

compressed  gas  whose  initial  parameters  were  P„  =  5  x  lO"*,  V_  =  0,  49  and  U  =0 

b  0  g 
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Fig.  1.  Pressure  versus 
Lagrangian  coordinate  R^. 


Fig.  3.  Density  versus  Rj. 


In  Fig.  4,  calculated  pressure,  velocity,  and  density  drop  at  the  shock  wave 
front  are  plotted  as  continuous,  dashed,  and  dot-dash  lines. 
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Fig,  4,  Shock  wave  parameters  In  water 
versus  time  t  £or  v  =  1,  2,  and  3, 


Simonov,  V,  A,  On  processes  originatinpi  during 
the  falling  of  an  impact  wave  on  a  tapered  cavity. 
FGiV,  no.  2,  1971,  280-284. 


Results  are  described  of  e}q>erimental  investigations  on  processes 
arising  from  an  impact  wave  of  subsonic  velocity  on  the  tip  of  a  tapered  cavity 
iu  metals,  and  a  study  »s  made  of  the  stability  of  the  flat  cumulation  jet  which 
is  produced.  E}q>eriments  were  conducted  with  aluminum,  copper  and  steel, 
the  arrangement  for  which  is  shown  in  Fig.  1.  The  cylindrical  explosive  charge 
of  diameter  40  cm  and  height  10-15  mm  consisted  of  type  TG  50/50  compound. 
Detonation  of  the  charge  was  synchronized  with  photographic  instruments  to 
obtain  a  true  picture  of  the  detonation  processes;  the  range  of  angle7(Fie.  1) 
was  varied  from  0.  5°  to  32®. 


Experiments  with  aluminum  specimens  showed  that  a  flat  cumulation 
metallic  jet  was  formed  during  the  normal  interaction  of  the  impact  wave  on  a 
blpered  cavity  with  an  angle  Zy  =  30O.  With  decrease  of  y.  jet  formation  dlaappears. 
that  is,  a  regime  of  collapsing  without  jet  formation  occurs,  characterized  by 
a  periodic  deformation  very  similar  to  wave  formation  during  explosion  welding. 

The  length  of  this  wave  was  found  to  vary  with  change  in  angle  y.  At  7=  1.  5°, 
wavelength  X  was  1.  5  mm;  at  7=  6.  50,  X  =  5.  5  mm.  In  the  range  7.6o.s.8o’  both 
wave  and  jet  regimes  were  found  to  exist  simultaneously.  At  7  ^  S®,  only  the 
jet  regime  existed,  and  at  7=150,  this  jet  was  capable  of  piercing  through  a 
duralumin  plate  of  10  mm  thickness. 


I 
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Fig.  1.  Experimental  sketch. 

I  -  metallic  specimen;  2  -  HI  charge 
with  flat  wave  regenerator;  3  -  detonator; 
4  -  metallic  plate  of  given  thickness. 


A  series  of  experiments  with  copper  samples  showed  a  very 
similar  process  to  that  with  Al,  and  the  transition  of  wave  Regime  to  jet  regime 
was  noted  between  the  same  interval  of  angles  7=  6°~8o.  The  wavelength  of 
aCu  sample  with  respect  to  7  is  given  below.  Increasing  7  from  2°  to  5°  increases 

T*  2  4  S  S.S  g  7 
1,  JTjf  t  2  2.5  3  3.5  3.5 

wavelength  from  1  to  3.  5  mm;  however  a  further  increase  in  angle  does  not 
increase  the  wavelength.  It  is  noted  that  geometrical  dissymmetry  does  not 
affect  the  angle  of  transition  from  wave  to  jet  formation,  and  that  the  jet  formed 
is  always  directed  along  the  bisector  of  angle  2  7  (Fig.  1). 

Similar  esqieriments  with  steel  samples  showed  that  transition  of 
wave  to  jet  regime  takes  place  at  an  interval  between  angles  7  =  40^50^ 

Experiments  were  also  conducted  with  samples  of  different 
materials,  but  having  a  like  geometry.  The  angln  of  transition  from  wave  to 
jet  regime  was  noted  to  increase  vory  significantly.  V'hen  one  cube  was  titanium 
and  the  other  copper,  the  jet  regime  was  noted  only  at  angles  7>20°;  and  the 
jet  was  inclined  from  the  symmetrical  axis  towards  this-  metal  having  a  heavier 
density,  in  this  case  copper. 

The  author  makes  a  comparison  of  his  results  with  those  of  Walsh 
et  al,  (Mekhanika.  Sbornik  perevodov,  vyp.  2(24),  1954),  but  they  do  not  seem 
to  be  in  good  agreement.  In  Walsh  et  al,  angles  of  transition  from  wave  to  jet 
formation  for  steel  and  Al  were  found  to  be  2  7  =  30°  and  2  7=  34°,  respectively. 
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taking  into  consideration  the  effects  of  incoming  flow  Uq.  In  ^e  present  article, 
U©  neglected^  as  it  depends  very  little  on  angle  7.  (U^  =  ^5,  y  ,  where 

D  -  velocity  of  impact  wave  in  sample,  U  -  body  velocity  of  particles).  According 
to  Walsh  et  al,  impact  of  solid  bodies  without  jet  formation  has  been  proven 
possible,  and  it  was  established  that  jet  formation  cannot  take  place  in  presence 
of  two  impact  waves  at  the  point  of  contact.  In  the  present  article,  the  author 
points  out  that  a  stationary  impact  wave  cannot  exist  at  the  point  of  contact  at  any 
angle  7,  as  the  falling  impact  wave  is  at  subsonic  speed,  so  that  the  jet  formation 
automatically  occurs.  However,  in  spite  of  the  above  discrepancies,  the  author 
still  agrees  with  the  possibility  of  the  process  occurring  without  jet  formation, 
but  with  generation  of  deformations  in  the  form  of  waves  on  the  collapsing  surface. 
He  points  out  that  in  such  a  case,  for  explaining  the  phenomenon  of  transition 
from  jet  to  wave  regime,  it  is  necessary  to  take  into  account  other  criteria 
related  to  the  model  of  an  ideal  fluid. 
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III;  SHOCK  WAVES  IN  GASES 


Pogorelov,  V,  I,  The  impart  of  an  incalculable 
supersonic  jet  on  a  plane.  IFZh,  v.  21,  no.  5, 

1971,  941-942. 

Shock  wave  separation  occurs  in  front  of  a  plane  positioned 
normally  to  the  nozzle  axis  when  an  incalculable  supersonic  jet  impinges  on 
the  plane.  The  separated  wave  is  convex  toward  the  plane  and  interacts  with 
tne  hanging  shock  Of  the  jet  if  the  nozzle-to-plane  spacing  is  such  that  the 
separated  wave  interferes  with  the  expansion  of  a  free  jet.  A  flow  measure¬ 
ment  equation  is  established  fot  gas  flow  ahead  of  and  behind  the  separated 
shock  wave  with  allowance  for  a  relationship  in  the  shock  wave  and  on 
assuniption  that  the  density  and  the  velocity  component  parallel  to  the  plane 
are  constant  behind  the  separated  wave.  This  flow  equation  is  used  to 
derive  a  differential  equation  which  defines  the  geometry  of  the  generatrix 
of  the  separated  wave.  The  boundary  condition  of  this  equation  is  the  derivative 
of  the  shock  wave  equation  for, a  triple  configuration  of  shock  waves.  The 
coordinate  of  the  triple  configuration  is  determined  by  a  simple  alg.braic 
correlation  using  the  numerical  value  of  the  derivative.  Gasdynamic  parameters 
of  subsonic  flow  behind  the  shock  wave  can  be  also  determined,  if  the  relative 
position  versus  the  plane  and  the  shape  of  the  shock  wave  are  known.  The 
calculated  coordinates  of  the  triple  point  and  geometry  of  the  separated  shock 
Y^ve  were  found  to  be  in  good  agreement  with  experimental  data  obtained 
from  flow  shadow  photographs. 


Pustovalov,  V.  K.  Self- similar  gas  motion 
behind  a  shock  wave  front  sustained  by  radiation. 

DAN  BelSSR,  v.  15,  no.  12,  1971,  1079  1081. 

Adiabatic  expansion  of  a  cold  perfect  gas  in  vacuum,  initiated 
by  a  strong  shock  wave,  is  analyzed  with  allowance  for  radiant  energy  q 
supplied  to  (q>0)  or  withdrawn  from  (q<0)  the  expanding  gas.  Shock  wave 
propagation  through  the  gas  is  sustained  by  radiation  which  determines  the 
motion  pattern  of  the  shock  front.  The  radiation  flux  density  is  defined  by 
q  =  qot*  ,  where  /  is  an  integer.  This  is  a  self- similar  gasdynamic  problem 
of  the  first  kind,  since  there  are  two  independent  dimensional  parameters 
q^  and  The  problem  is  described  by  the  set  of  equations 


(1) 


f"' it""’ 

+  fc.-a- 


where  Pi,  pi,  vj  are  the  unknown  dimensionless  density,  pressure,  and 
velocity  orthe  gas,  f  is  a  self- similar  variable,  Y  is  the  adiabatic  exponent, 
«  =  1  +  //3  IS  the  self- similarity  index,  k  =  f/S,  and  m  =  2  //3.  The 

oun  ary  conditions  at  the  shock  frontand  on  the  vacuum  side  complete  the 
formulation  of  the  problem. 

Solution  of  (1)  with  the  boundary  conditions  is  given  as  follows: 
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where  a  is  one  of  the  dimensionless  constants  which  define  the  boundary 
condition  at  the  shock  front.  Solution  (2)  for  a  =  0  and  Y=  7/5  describes 
shock  wave  propagation  in  the  absence  of  radiation  (q  =  0).  Solution  (2) 
lor  a  >0,  i.e. ,  q>  0,  describes  propagation  of  a  shock  wave  with  radiant 
energy  absorption.  A  numerical  solution  of  (.1)  with  boundary  conditions 

lor  a<0,  i.e. ,  q<0,  describes  propagation  of  a  shock  wave  front  with 
self- radiation. 


Afanasenkov,  A.  N. ,  I.  M.  Voskoboynikov,  and 
A,  Ya,  Apln.  Transmission  of  detonation  through 
an  air  gap.  IN;  Detonatsiya  vzryvchatykh  veshchestv 
i  bezopasnost'  vzryvnykh  rabot.  Moskva.  Izd-vo 
Nedra,  1967,  lOl-llO. 

The  purpose  of  this  work  was  to  produce  a  simple  criterion 
for  explosive  sensitivity  and  to  study  the  possible  mechanism  of  transmission 

of  detonation  through  an  air  gap  on  the  basis  of  available  data.  Experiments 
were  performed  using  passive  charges  20  mm  in  diameter  and  70  mm 
long  of  hexogene  of  various-densities,  and  active  charges  of  hexogene  and 
TG  50/50  of  various  weights  and  shapes.  The  experiments  determined  the 
significance  of  explosion  products  and  the  air  shock  wave,  and  produced  a 
quantitative  characteristic  of  the  sensitivity  of  explosives  to  detonation 
transmission  through  an  air  gap.  The  dependence  of  the  distance  of  trans¬ 
mission  of  detonation  on  explosive  charge  density  was  explained. 
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IV»  SHOCK  WAVES  IN  LIQUIDS 


Kozachenko,  L.  S.  and  B.  D.  Khristoforov. 
Parameters  of  a  shock  wave  from  explosion 

17?!^  ^  FGiV,  no.  1 
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the  experimental  condillons  Th«I  •'  7  eliminated  under 

in  separate  vessels  with  ,"r  f  '*Periments  were  carried  out 

concrete  bottoms  (ill)  The  three"^bott  P°ly»tyrene  (11),  or 

Ot  longitudinal  and  C  ot  transvers^^^^^^^ 

types.  Depth‘charBes^;f  100  ka  c^ift -Nt  T'*""' 
compacted  TNT  were  detonated  in  3  1  L  f  \ZT 
bottoms,  respectively,  at  relative  distances  H  =  1.11  uu  .  1 
radu  from  the  bottom.  The  pressure  sensors  were  locLd  at  fd!mh  b 

for  vess'enl  c^rgT"' 

21  anH  fK  oscilloscope  traces  of  shock  waves  in  water  (Fie  1  anH 

2)  and  the  experimental  T(h)  and  t(  e  ,  plots  (Fig.  3  and  4)  are  shown  The 

rnterl^cr"^^°!h‘'  ^^gg^sts  a  nonlinear 

decreasrthrf^^^^  Propagating  from  the  bottom.  The  latter 

decrease  the  frontal  pressure,  as  shown  by  the  plots  of  p(h).  In  the  case 

obse^  ^  type  bottom,  p^^^x  of  the  reflected  compression  wave,  which  is 

attainr30%"of  ^in^^he^  1^"^^^  of  incidence  a  of  the  forward  wave  to  the  bottom, 
attains  JO/o  of  p  m  the  leading  shock  wave.  A  bottom  wave  behind  and  a  learlin 

shock  wave  ahead  of  the  forward  wavefront  were  observed  in  the  case  of  the 
II  type  bottom  (Fig.  2).  The  leading  wave  is  generated  only  at  an  angle  «  ~  « 
which  IS  determined  by  the  intersection  of  two  wavefronts.  The  p  ^ 

exceed  by  40%  the  explosion  pressure  in  an  infinite  liquid  volume  die^'to  the 
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Fig.  1.  Relative  pressure  p/pj  versus 
reduced  duration  T^sec/m)  of  explosive 
shock  waves  in  a  water  basin  with  a  sand 
bottom. 


Fig.  2.  Relative  pressure  p/pj  versus 
reduced  duration  T  (msec/m)  of  explosive 
shock  waveL>  near  the  polystyrene  bottom. 


Fig.  3.  Duration  t  (msec/m)  of  shock 
wave  versus  distance  h  from  the  pickup 
to  the  sand  bottom. 


Fig.  4,  Reduced  duration  T  (msec /m) 
of  a  shock  wave  versus  angle  «  for  an 
explosion  at  a  polystyrene  bottom. 


superposition  of  the  leading  and  forward  wave  pressures  (Fig.  Z,  d,e).  The 
Pmax  •  “  •o*  20-30°  increases  sharply  because  of  the  interaction 
between  the  forward  and  leading  waves,  as  shown  on  the  p(*)  plot,  Interacticn 
of  the  forward  wave  with  the  bottom-reflected  rarefaction  wave  manifests 
itself  in  a  decrease  of  t  with  a  decrease  in  «  (Fig.  4).  The  shaded  area  in 
Fig.  4  designates  the  near-bottom  zone,  where  total  t  of  the  wave  field  is  much 
longer  than  t  of  the  forward  wave,  i 

Tabulated  experimental  data  (Table  1)  show  that  the  parameters 


/?=90,  //=•! 

/?b90,  //=.3 

0 

0.28 

0.45 

0.384 

0 

0.38 

0.40 

1 

0.21 

0.64 

O.S-H 

0.123 

0,70 

0.55 

2 

0.78 

0.75 

0.376 

0.25 

0.78 

0.61 

3 

0,73 

0.78 

0.376 

0.37 

.0.99 

0.72 

7?=60.  W=1 

/?=120.//=1 

0 

0.39 

0.31 

3.5 

0 

0.33 

0.30 

4 

0.52 

0.58 

2.1 

0  49 

0.41 

0..51 

tt 

0,57 

0.71 

0,93 

0.87 

0.46 

0.62 

1 


Table  1.  Parameters  p/pj  and  t  of  explosive 
shock  wave  in  a  water  basin  with  concrete 
bottom.  Footnote:  The  indexes  e  and  p  indicate 
the  experimental  data  and  the  data  calculated 
from  the  formulas  given  in  the  text,  respectively. 


23 


of  the  forward  wave  for  the  III  type  bottom  are  greatly  Influenced  by  the 
bottom-reflected  rarefaction  wave.  The  p^ax  forward  wave  decreases, 

when  the  sensor  is  positioned  closer  to  the  bottom,  and  approaches  the  Pmax 
value  calculated  using  an  empirical  formula.  However,  the  experimental  ^ 

T  disagrees  with  the  t  calculated  from  an  empirical  formula,  because  of  the 
effect  of  a  strong  bottom  wave. 

The  experimental  data  show  that,  at  great  distances  from  the 
blast  center,  the  effect  on  shock  wave  parameters  in  water  of  all  three  bottom 
types  is  analogous  to  that  of  a  free  surface.  The  empirical  formulas  proposed 
can  be  used  to  calculate  shock  wave  parameters.  The  actual  t  of  the  shock 
wave  measured  near  the  bottom  is  usually  longer  than  the  calculated  t  ,  because 
of  bottom  wave  interference.  Pressure  in  the  lateral  and  leading  waves  is 
generally  not  over  20-30%  of  that  in  the  forward  wave.  The  maximum  pressure 
near  the  intersection  of  the  leading  and  forward  waves  is  30-40% 
higher  than  the  explosion  pressure  measured  in  an  infinite  liquid  volume. 

This  paper  is  an  extension  of  a  similar  work  by  the  authors  reported  in 
Effects  of  Strong  Explosions,  no.  2,  p.  111. 

Podurets,  M.  A.,  G.  V.  Sumakov,  R.  F.  Trunin, 

Popov  and  B.  N.  Moiseyev.  Compression  of 

water  by  strong  shock  waves.  ZhETF\  v.  62,  no.  2. 

1972,  710-712. 

Measurement  of  shock  compressibility  of  water  by  absolute 
methods  is  limited  to  a  maximum  pressure  of  1~2  Mbar.  For  higher  pressures, 
compressibility  can  be  found  only  by  relative  methods,  based  on  a  reliable 
interpolation  of  the  impact  adiabatic  curve  of  a  standard  analog  substance,  which 
in  case  of  water  can  be  aluminum.  In  this  article  the  authors  use  the  cited  method 
to  investigate  the  relative  compressibility  of  water  at  pressures  ~14  Mbar.  The 
measurements  were  obtained  bv  successive  recording  of  a  shock  wave  propagating 
through  an  aluminum  block  of  thickness  A=  160  mm  and  a  water  layer  of  thickness 
A  =  80  mm;  wave  velocity  (D)  was  determined  with  an  error  <  1%.  Parameters 
of  compressibility,  according  to  this  known  wave  velocity,  were  determined  by 
P-U  diagrams  (pressure  -  body  velocity  of  the  substance  behind  shock  wave  front). 
Damping  of  the  shock  wave  at  the  interface  between  aluminum  and  water  surfaces 
is  taken  into  account  by  a  small  correction  in  the  calculations.  Initial  parameters 
in  aluminum  are  found  to  be  =  36.  40  km/sec;  Uai  =  25.  55  km/sec;  Pai  =  25.20 
Mbar,  and  for  water  they  are  D  =  43.95  km/sec,  U  =  32.42  km/sec,  P  =  14.25 
Mbar,  and  p=  2.815  g/cm  .  Errors  in  determining  the  water  density  are  connected 
with  the  experimental  inaccuracy  in  measuring  the  wave  velocity  in  both  media, 
and  also  with  the  uncertainty  in  the  standard  adiabatic  curve  for  aluminum;  they 
amount  to  about  ^e/eteO.  03.  Results  are  plotted  in  Fig.  1.  The  figure  includes  a 
curve  calculated  according  to  the  quantum- statistical  Thomas-Fermi  model  A 
comparison  of  calculated  and  extrapolated  experimental  dynamic  adiabatic  curves 
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Fig.  1.  Shock  adiabati  of  woter. 

A  -  author  a*  reaulta 

I  -  experimental  adiabats,  other  authors. 

of  water  shows  the  applicability  of  the  Thomas-Fermi  model  for  determining 
the  parameters  of  compressibility  for  relatively  light  substances  such  as  water 
In  the  region  of  very  high  pressures,  P^400  Mbar. 


Ocheretin,  V.  N.  Study  of  shock  wave  formation 
from  explodin;;  wires  in  water.  EOM,  no.  5 
1971,  39-41.  ' 


I,  ,  ^  ‘*®®c>’lptlon  is  given  of  an  experiment  on  the  formation  of 

iven  .:fu  .vetch  U 

Ji  j  ^  conducted  with  steel  and  copper  wires  of  0  15  mm 

Mpaclty  of  1/tf.  Photograph,  war.  taken  at  the  rate  of  600,000  frame, 
per  second;  two  photographs  are  given  illustrating  steel  and  copper  wire 
..pl..lon..  It  wa.  noted  that  the  channel  formed  from  the  expIoTon  Ta 
Steel  ^re  had  a  uniform  structure,  while  that  formed  by  copper  wire  was 
Energy  levels  in  the  channels  were  measured  at  150  kj/cm^ 

****^  copper  wire  respectively.  From  the  test 
results,  the  author  recommends  use  of  steel  wires,  rather  than  copper  or 
brass,  for  its  superior  hydrodynamic  parameters. 
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high  voltage  step-up  transformer;  V  - 
rectifier:  IVPl,  IVP2  -  spark-gaps; 

B  •  test  vessel;  ♦  -  aluminum  foil; 

BP  -  exploding  wire. 
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y,  SHOCK  WAVES  IN  SOLIDS 


Limarev,  A,  Ye,  and  A,  D,  Chernyshov. 

Propagation  of  shock  waves  in  a  reinforced 
elasto- plastic  medium.  PMM,  no.  6,  1971, 

1083.1088. 

Propagation  and  structure  of  shock  waves  in  a  translationally 
-reinforced  elasto-plastic  medium  is  described  by  the  closed  set  of  discon¬ 
tinuity  equations 


Io,/*I  -  2|i,  =  2  I<pl  {2n,n,c*,^  _  2|i,3*4j  -  p,p,  + 


I®»kl  “  (31.1  +  2fi|)  (2) 

[»»•  -  vj{2  (.„•  -  .„•)•  -  [v  -  »,„.]!  -  0  ,3, 

and 

V/  =.  —  p*G  (v,l,  (u,,  /I  .  (o,v/,  (»>,)  or  —  Cw,  (4) 

where  o^j,  o  are  stress  coniponents,  e^j  are  strain  components, 

/^o'  /^l'  /*Z  elastic  coefficients,  v  j  is  normal  to  discontinuity  surface, 

G  is  the  propagation  velocity  of  a  nonremovable  discontinuity  surface  in  a  , 
continuous  medium,  Uj,  Vj  are  the  displacement  and  velocity  of  a  particle, 
and  the  asterisk  denotes  a  tensor  deviator  component.  Equations  1-4  were 
derived  assuming  a  linear  dependence  between  two  discontinuous  functions 
within  the  shock  wave  transition  layer.  This  linear  approximation  made 
possible  a  common  solution  to  the  problems  of  propagation  and  structure  of 
shock  waves  in  a  plastic  noncompressible  raaterial. 

Using  (1)  and  (4),  an  expression  is  derived  for  the  characteristic 
of  irreversible  plastic  deformation  of  a  continuous  medium 

2  (7) ».  f,(o„ 

IFt  -  p*c«  -  Xt  -;2p„  F,  -  p*G» X,  _ »/,  p,  (1 .}.  ,,,/p,) 

where  is  the  stress  component  of  a  plastic  element  p  of  the  rheological 
model  of  the  material,  the  index  n  and  the  sign  +  indicate  direction  of  the 
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normal  v  ^  and  the  value  p  ahead  of  the  shock  wavefront,  respectively.  The 
velocity  p'^G^  of  discontinuity  surface  propagation  is  delerr.iined  by  the  roots 
of  the  equation  of  the  4th  degree  in  p'*'G2,  which  was  derived  from  the 
equations  1-5.  The  root  p'^G^  cited  equation  corresponds  to  a 

transverse  neutral  wave  without  discontinuity  of  plastic  deformation  and 
satisfies  the  condition  f  J  =  0.  This  particular  case  is  illustrated  by  an 
example  of  a  rectangular  elasto- plastic  plate  subjected  to  either  a  tensile  or 
a  shear  stress.  Two  other  roots  of  the  4th  degree  equation  indicated  the 
possibility  of  generating  two  plastic  shock  waves  in  a  three-dimensional  re¬ 
inforced  elasto-plastic  medium.  One  of  the  plastic  waves  propagates  in  a 
medium  under  spherically  symmetric,  stress,  the  other  in  a  medium  where 
S|«  =  0. 


Godunov,  S.  K. ,  A.  A.  Deribas,  I.  D.  Zakharenko, 
and  V.  I.  Mali.  Investigation  of  viscosity  of  metals 
during  high  velocity  shock.  FGiV,  no.  1,  1971,  135-141. 

This  article  describes  a  study  of  metal  viscosity  during  shock 
from  explosive  welding.  An  experiment  is  described  in  which  steel  plates  were 
welded  as  shown  in  Fig.  1.  On  detonation  of  explosive  (5),  the  upper  plate  (3) 
is  driven  onto  the  lower  plate  (1),  as  a  result  of  which  the  two  plates  are  welded 
as  shown  in  Fig.  lb.  The  welded  specimen  was  then  cut  along  the  direction  of 


Fig.  1.  Explosive  weld. 

1, 3  -  steel  plates;  2  -  slot;  4  -  hole 

(0.  3~0.  5  mm);  5  -  explosive;  6  -  detonator; 

7  -  wooden  base. 

movement  of  the  contact  point  an^  the  cross-section  was  microphotographed 
(Fig.  2).  Horizontal  displacement  z  with  relation  to  the  distance  y  from  the 
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Interlace  wai  meaiured  on  the  photograph.  Aty>4j(<j.  thlckneii  of  upper 
plate,  <2  -  thlckneii  of  lower  plate),  the  relatlonihip  between  *  and  y  followi  a 
parabolic  curve,  z  =  a(y  -  «2)2;  while  at  y  <  *j,  z  =  -  an  exponential 

curve  (a,  b  =  cunitanti).  Dliplacement  of  Impact  polnti  of  the  platei  with 
relation  to  metal  vlacoslty  wai  Inveitlgated  and  waa  found  to  be  Inver aely  p^o*. 
portional  to  viacoilty.  Mathematical  expreaalona  are  derived  for  determining 
■tored  Impultiei  and  the  diffusion  of  horizontal  velocltiea  In  both  upper  and  lower 
platei.  Using  experimental  data  and  derived  formulae,  vlecoelty  coefflclente 
were  determined  for  Al,  Cu  and  iteel  (Table  1).  Vlecoelty  coefficient!  for  Al, 

Cu  and  itesl  are  leen  to  have  an  Increasing  order  of  (0.  3M).  8)lo5p,  (2~  2.  7)105p 
and  (4'-5)l{  5p  respectively. 


Fig.  2.  Microphotograph  of  weld  section. 

1  -  lower  plate;  2  -  compressed  plate; 

3  >  upper  plate;  4  -  compressed  wire. 

Arrow  Indicates  the  direction  of  detonation. 

The  authors  compare  these  values  with  those  of  Popov  (Inzhenernyy 
sbornlk,  1941,  1,  1)  and  Il'yuihln  (Inzhenernyy  sbornlk,  1941,  1,  1),  and  the 
value  of  viscosity  coefficient  for  steel  Is  found  to  bo  the  same  In  both  cases. 
Similarly,  the  obtained  value  for  Al  is  found  to  be  In  good  agreement  with  the 
results  of  Sakharov,  Zaydel',  et  al  (DAN  SSSR,  1964,  159)  and  Mlneyev  and 
Savlnov  (ZhETF,  1967,  57,  3).  However,  according  to  the  results  of  Il'yuihin 
(Uchenyy  zap.  MGU,  Mekh,  1940,  39,  II  and  Inzhenernyy  sbornlk,  1941,  1,  1), 
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aluminum  and  steal  have  approximate!/  the  same  coefficient  of  viscosity  (about 
lO^p),  which  lies  far  below  the  results  of  Table  1.  Similarly,  works  of  Sakharov, 
Zaydel',  etc,  as  well  as  Mineyev  and  Savinov,  disagree  with  the  present  results 
in  the  fact  that  in  their  results  Al,  Cu  and  steel  all  have  the  same  viscosity 
value.  The  authors  discuss  these  discrepancies,  and  express  doubts  as  to  the 
accuracy  of  the  cited  other  experiments. 


Materials 

ll. 

tM 

€M 

u 

M  e 

T 

►•ID”*-  poise 

]r-4L4 

jreOJ 

Al  (  n  16) 

Cu  (M  3) 

Steel^^’’ 

0.4 

0.4 

0.4 

0.4 

0.45 

0.4 

2.0 

3.0 

1.4 

2.4 
2.8 
3.0 

2.5 

3.1 
1.7 

4.2 
3.1 
4.0 

10*30' 

14* 

21* 

14*20' 

15*20' 

14*20' 

0.31 

0.81 

2.1 

2.5 

3.9 

0.32 

0.86 

2.7 

2.0 

4.1 

4.8 

Table.  1.  Shock  viscosity  results. 


Zubarev,  V.  N.  Shock  wave  compression  of  piezoelectric 
ceramics.  ZhPMTF,  no.  2,  1971.  106-110. 


Dynamic  compressibility  and  piezoelectric  characteristics  of  type 
TsTS  53/47  lead  zirronate  titanate  ceramic  samples  were  measured  at 
pressures  in  the  100-500  kbar  range  in  front  of  a  shock  wave.  Also,  nonre- 
versible  changes  wore  obaerved  in  the  samples  after  their  short-time 
compression  at  350  and  500  kbar.  Shock  waves  from  detonation  of  an  explosive 
charge  were  transmitted  to  the  sample  through  a  metallic  screen.  The 
measured  wave  velocities  U  in  the  sample  and  the  calculated  mass  velocities 
u,  pressures  p,  and  temperatures  T  were  tabulated.  An  experimental  plot 
of  p  versus  specific  volume  V  shows,  in  contrast  to  earlier  American  data, 
the  absence  of  a  two-wave  shock  structure  in  the  100-350  kbar  range;  tne 
linearity  of  the  U  -  u  dependence  confirmed  this  finding.  The  measurement 
of  current  during  shock  wave  transmission  Indicated  only  small  variations  in 
polarization  at  a  load  ^  1  ohm,  so  that  inductance  Dj  can  be  approximately 
described  by  Dj  =  +  CqE.  It  follows  that  dielectric  constant  c  and  electric 

conductivity  X  in  a  compressed  sample  can  be  calculated  from  the  equation 
of  R.  M.  Zaydel'  [ZhETF,  v.  54,  no.  4,  1968]  using  the  experimental  extreme^ 
values  of  current.  Calculations  revealed  that  X  is  relatively  high  4.  5  x  10" 
ohm"^  X  cm-1)  at  p  =  105  kbar,  and  increases  by  a  fac  -^or  of  about  20  at 


30 


°«n,T3Vto 

to  The‘p.'r«l^Src  *‘h^'  f  “''"‘t  •'“  “>  ‘''«  conclu.’ior[hlt  r.ion 

rano*  A  1  phase  is  possible  at  shock  wave  pressures  in  the  300  -  400  kbar 
range.  An  analogous  experiment  by  Zubarev  et  al  on  shock  compression  of 
plastics  was  reported  in  Explosion  Effects  Report  no.  2,  p.  4. 
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VI.  FLOW  PAST  BODIES 


Balakin,  V.  B.  and  V,  V,  Bulanov.  A _ 

numerical  solution  of  the  problem  of  interaction 
of  a  shoe’  wave  with  a  cylindrical  body  in 
supersonic  flow.  IFZh,  v.  32,  no.  6,  1971, 

1033-1039. 

A  two-dimensional  finite  difference  scheme  of  the  second  order 
accuracy  is  introduced  to  compute  the  process  of  collision  of  a  nonstationary 
shock  wave  with  a  stationary  flow  around  the  nose-end  of  a  cylinder  moving 
at  a  supersonic  speed.  The  initial  set  of  equations  of  two-dimensional  flow 
in  divergent  form  is  presented  in  vector  notation 

+  —ll{l)-‘0  (1) 

di  Ox  dtj  y  '  ’ 

and  allowance  is  made  for  through  computation  of  shock  waves.  Eq.  (1)  is 
then  approximated  in  two  steps  by  a  finite  difference  equation  which 
contains  artificial  viscosity  terms  in  the  x  and  y  directions  measured 
along  and  from  the  flow  axis,  respectively.  The  necessary  and  sufficient 
stability  criteria  of  the  finite  difference  scheme 


were  established  by  a  Fourier  method.  A  =  dF/df  and  B  -  dG/df  in  (2) 
are  Jacobians  of  the  F  and  G  vectors  in  a  random  point  of  the  flow, A  *  and 
Ay  are  the  difference  grid  pitches  along  x  and  y,  A  t  is  the  time  pitch,  and 
E  is  the  matrix  unit.  By  use  of  a  through  computation  method,  a  commong 
numerical  solution  is  outlined  for  two  corollary  problems?  the  initial 
stipersonic  streamline  flow  around  the  body,  and  the  oncoming  shock  wave 
with  establishment  of  a  new  regime  of  supersonic  streamline  flow.  The 
initial  flow  parameters  are  those  of  the  first  shock  wave.  The  parameters  of 
the  second  shock  wave  are  calculated  by  means  of  Hugoniot  formulas.  Numer 
ical  solutions  of  the  problem  are  illustrated  graphically. 


Poluboyarinov,  A.  K.  Propaeation  of  a 
shock  front  reflected  by  a  blunt  body. 

MZhiG,  no.  2,  1971,  70-77. 

Shook  wave  front  interaction  v'ith  a  blunt  sphere  or  cylinder  is 
analyzed  on  the  assumption  thS't  the  front  plane  is  normal  to  the  symmetry 
axis  or  plane  of  the  body  and,  in  the  case  of  the  cylinder,  is  parallel  to  the 
body  generatrix.  An  exact  differential  equation  is  derived  to  describe 
nonstationary  separation  of  the  symmetry  center  of  the  shock  reflected  by  a 
static  body.  The  equation  describes  dimenlsionless  axial  shock  separation  c 
as  a  function  of  the  independent  variable  M  of  the  reflected  front.  The  equation 
also  contains  two  additional  unknown  quantities:  the  derivative  {  dv/  d  x)i  of 
gas  velocity  immediately  behind  the  shock,  and  the  radius  r  of  the  reflected 
shock  front  curvature  along  the  symmetry  axis.  An  exact  formula  was 
derived  for  initial  (  d  v/  3  S'S  a  function  of  r^  and  Approximate 

formulas  define  r,j(  and  (dr/dT),^.  in  the  initial  separation  phase  of  the  front 
from  the  body.  A  fairly  simple  approximate  method  is  introduced  to 
determine  nonstationary  e  for  a  sphere,  on  the  basis  of  the  derived  formulas 
and  ♦'he  precise  equation  of  £ .  The  method  consists  of  integratinr,  the 
equation  d  £ /dM  =  V>,  where  'Pis  a  function  of  r  and  (av/3x)j.  A  numerical 
example  shows  that  the  method  correctly  describes  the  initial  separation 
phase,  and  also  its  final  phase  if  a  correction  to  the  ^function  is  introduced. 

Yeliseyev,  Yu.  B.  and  A.  Ya.  Cherkez.  The 
effect  of  an  increase  in  stagnation  temperature 
in  a  streamline  gas  flow  past  a  deep  cavity. 

M Z hiG,  no.  3,  1971,  8-18.  ^ ' 


fKo  *  earlier  observed  anomalous  increase  in  temperature  of 

the  separated  gas  flow  within  a  cavity  above  the  stagnation  value  To  is 
iscussed,  on  the  basis  of  a  theoretical  model  of  compression  and  decom¬ 
pression  waves  in  a  cylindrical  tube.  Piston  oscillatory  movements  at  the 
open  end  of  the  tube  model  simulate  conditions  at  the  streamline  flow- separated 
gas  interface.  Characteristics  of  the  nonstationary  self-os^ilTat^^^^^^^^^ 

givin  ^To^r  *  interface  are  formulated  and  a  numerical  example  is 
g  en.  A  preliminary  conclusion  was  drawn  that  the  observed  phenomenon 

Incre^sed^r  dimensions  but  is  influenced  by 

models  at  M°<"l  and  M  "  experiments  with  different 

..  *4  1,  ^  ^  M  -  1.6  and  2.0  are  described.  In  agreement  with 

modef  calculations,  the  stationary  air  tempereture  at  the  base  of  a 

“  4irmerrh minimum  of  100°K  above  T„.  The  experimental  data 
CO  firmed  the  .hypothesis  that  the  shock  wave  thermal  effect  and  heat  release 
ollows  energy  dissipation  in  shock  waves.  The  intensity  of  the  observed 


bv  ^*^8®  and  the  angle  formed 

possibility  of  reproducing  the  effects  in  the  surface  cavity  of  an  obfe^ct  flvine 
are  gi\  the  undesirable  thermal  effect  ^ 


Bezhanov,  K.  A.  Diffraction  of  a  shock 
wave  on  a  wedge  moving  at  supersonic 
velocity.  Prikladnaya  matematika  i 
mekhanika,  no.  4,  1969,  631-637. 

.tr„n.  V  ■*  ““'*y  "«<<«  o£  the  diffraction  of  an  arbitrarily 

Ualne  th^  "'“'’‘"S  ““P'^onic  velocity 

ng  the  notations  of  Fig.  1,  the  author  assumeb  that  the  incident  angle  8of'^ 


Fig.  1.  Flow  model. 

no'tnH  "  the  as  the  wedge  angle  a.  In  this 

notation  Uq  defines  absolute  shock  wave  velocity,  while  U  is  its  velocity 

8o»"e  simplifying  assumptions 

the  analysis  is  made  with  a  linearized  equation  for  planar  gas  flow,  uf.ing 
Imensionless  parameters.  The  work  is  an  extension  of  some  cited  earlier 

treatments  of  a  similar  problem  in  which  diffraction  of  the  shock  wave  was 
assumed  not  to  take  place. 


Borovoy,  V.  Ya.  and  M.  V.  Ryzhkova. 

Heat  transfer  on  a  semi-conic  surface  at 
high  hypersonic  velocities.  MZhiG.  nn,  d 
1969,  137-142. 

Formation  of  narrow  bands  of  intensive  heat  transfer  on  the 
^  streamlined  spherically-blunted  semi-cone 
7  A  in6®  ®*P«'i^entally  at  a  Reynolds  number  Re  in  the  1.  5  x  105 

-  2.  4  X  10  range.  At  a  sufficiently  high  Re,  the  boundary  layer  separates 


Fig,  1.  Spherically  blunted  semi-cone. 

from  the  flat  surface  and  two  separate  vortex  formations  appear  near  the 
edges  as  boundary  lines  along  the  surface.  Most  of  the  heating  occurs 
along  the  boundary  lines  where  heat  transfer  coefficient  h°  may  be  several 
times  higher  than  h°  on  the  remaining  part  of  the  flat  surface.  The 
experiments  show  that  Re  is  the  most  important  factor  in  formation  of  the 
boundary  lines  of  intensive  heat  transfer.  At  Re  =  2.  5  x  lo4  -  10  x  10^, 
^max  increases  nearly  proportionally  to  Re,  while  at  Re  >  10  x  10"^  the  rate 
i'max  increase  is  somewhat  slower.  The  strong  effect  of  Re  confirmed 
the  assumption  that  thelaminar  flow  changes  to  turbulent  at  the  boundary  of 
the  separation  zone.  Boundary  layer  separation  occurs  because  of  an  over¬ 
flow  of  gas  from  the  conic  to  the  flat  surface.  This  was  confirmed  by  an 
experiment  with  a  barrier.  Experiments  with  a  reduced  thickness  boundary 
layer  and  increased  conic  surface  roughness  did  not  contradict  the  proposed 
flow  pattern. 


Myshenkov,  V.  I.  Subsonic  and  tfan sonic 
viscous  gas  flow  in  the  wake  of  a  two-dimensional 
body.  MZhiG,  no.  2,  1970,  73-79. 

Viscous  gas  flow  in  the  wake  of  a  finite  thickness  plate  is 
analyzed  at  subsonic  and  transonic  velocities,  Prandtl  number  =  0.71,  and 
Reynolds  numbers  R  in  the  1-1,000  range.  The  problem  is  formulated  by 
Nayier-Stokes  equations  for  given  flow  parameters  at  infinity  and  at  a 
finite  distance  from  the  bottom  of  the  plate.  Equations  of  state  and  motion 
are  solved  by  approximation  of  a  finite-difference  scheme  in  two  steps. 

A  nurnerical  solution  which  satisfies  the  condition  of  a  steady  flow  is  presented 
graphically  for  different  R  and  M  values.  The  solution  shows  that  at  R<  1.  7 
gas  flows  without  separation  from  the  bottom,  and  at  R>1.  7  a  region  of  reverse 
flow  with  separation  forms  at  the  stagnation  point.  The  reverse  flow  region 
expands  with  increase  in  R  and  the  separation  point  shifts  upwards  in  the 
direction  of  the  angular  point  of  the  plate.  A  deformation  of  the  reverse 
flow  region  with  a  decrease  in  relative  velocity  of  reverse  current  is 
observed  at  R  =  100  at  increased  M  (from  0.288  to  1. 15).  The  flow  patterns 
are  in  agreement  with  earlier  theoretical  and  experimental  data.  Pressure 
variation  patterns  at  the  bottom  of  the  plate  are  also  shown  in  the  flow 
plane  and  the  plane  normal  to  flow. 
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Stulov,  V.  P. ,  and  L.  I.  Turchak.  Non¬ 
equilibrium  chemical  reactions  in  a  shock 
layer  in  a  streamline  flow  of  a  carbon  dioxide 
-nitrogen-argon  mixture  around  a  sphere. 

MZhiG,  no.  5,  1969,  147-150. 

A  theoretical  evaluation  of  hypersonic  (M  =  16-20)  flow  of 
C02-N2-Ar  mixtures  around  a  sphere  was  made  using  a  known  physico-chemical 
model  of  a  shock-heated  gas  mixture.  The  model  presumes  the  existence 
of  nine  components  in  the  shock  layer,  with  allowance  for  five  molecular 
dissociation  and  five  atom  exchange  reactions  which  are  tabulated.  Dissociation 
rate  constants  were  calculated  for  collisions  with  particles  of  different 
catalytic  activity.  The  presence  of  O2  was  disregarded  in  the  establishment  of 
the  model,  and  molecular  vibrations  at  T  S  4000°K  were  considered  to  be  at 
equilibrium.  It  was  shown  that  the  flow  changes  significantly  with  a  change 
in  relative  concentrations  Yoo  (CO2)  and  (N2)  in  the  shock  layer.  Shock 
wave  separation  increases  by~50%  when  Y^o  (CO2)  is  decreased  from  0.95 
to  0.05,  and  thermal  energy  absorbed  by  the  gals  simultaneously  increases 
because  of  a  decrease  in  T  and  an  increase  in  the  fraction  of  atomic  components. 
The  effect  of  Ar  addition  is  analogous  to  that  of  N2.  Most  of  the  heat  is 
absorbed  by  CO2  dissociation  because  of  its  higher  rate. 


Shifrin,  E.  G.  Formation  of  a  ’’hanging" 
compression  shock  in  the  streamline  flow 
past  a  profile  with  a  broken  generator.  PMM, 
no.  6,  1970,  1159-1167. 


Uniform  supersonic  streamline  flow  past  a  profile  with  a 
convex  salient  point  is  analyzed  in  the  case  of  an  oncoming  stream  with  a 
low  supersonic  velocity,  when  entropy  changes  in  the  shock  wave  can  be 
disregarded.  The  region  behind  the  shock  wave  is  projected  on  the  hodograph 


Fig.  1.  Projection  of  the  flow  behind 
shock  wave  on  hodograph  plane  T]  -  /3  . 
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plane  tj  (Fig.  1)  where  T|=  (k  +  l)l/3(  p.  .  i),  ^  the  inclination  of  the 
velocity  vector  to  the  symmetry  axis,  k  is  the  velocity  factor,  and  k  is 
the  adiabatic  exponent.  A  streamline  flow  with  attached  shock  wave  past 
a  profile  formed  by  straight  line  intercepts  in  considered  in  a  first  approxi¬ 
mation.  The  equations  of  the  shock  polar /3  and  intersecting  characteristics, 
e.  g. ,  aj  are  given  in  a  transonic  approximation.  Proof  is  offered  that 
a  single- "hanging"  shock  is  necessarily  generated  behind  the  shock  wave. 

The  "hanging"  shock  propagates  downstream  from  its  starting  or  end  point 
which  is  located  at  an  infinite  distance  from  the  profile.  The  shock  at  /3  ^  0 
is  convex  in  relation  to  the  flow  behind  it.  In  addition  to  the  tail  shock,  there 
is  an  internal  shock,  at  /9o<0,  which  intersects  the  tail  shock  and  forms 
a  lambda. shaped  shock. 


Antonets,  A.  V.  Computation  of  three 
dimensional  supersonic  flow  past  blunt  bodifta 
with  generator  breaks  and  allowance  for 
equilibrium  and  frozen  st^es  of  the  gas  in  the 
shock  layer.  MZhiG.  nn.  7.^  1  Q7n  170 

sonir  /  An  algorithm  is  developed  for  computation  of  the  steady  super¬ 

sonic  flow  of  an  inviscid  gas  past  a  blunt  body  with  multiple  breaks  of  a 
generator.  The  algorithm  is  based  on  the  principle  of  separation  (behind 

each  of  an  arbitrary  number  of  breaks)  of  the  regions  of  smooth  solutions 
between  the  body  and  the  eharacteristic  shock  wave  surfaces,  at  a  positive 
(expansion)  or  negative  (compression)  break.  Finite-difference  equations 
In  the  nodal  points  of  a  net  are  substituted  for  differential  equations  in  the 
regions  of  smooth  solution.  The  unknown  velocity  components,  pressure  p, 
and  enthalpy  i  are  determined  from  the  equations  by  the  method  of  successive 
calculations.  The  algorithm  accounts  for  the  effect  of  nonequilibrium 
physicochemical  reactions  by  computing  the  flow  for  either  the  entire  sho-k 
layer  in  a  nonequilibrium  (frozen)  state  or  for  only  a  partially  frozen  layer 
in  the  area  adjacent  to  the  body  surface,  the  other  part  being  occupied  by  a 

fs‘ ilLslr^tld  I  ;  i  The  effect  of  the  gas  equation  of  state  in  the  shock  layer 
aoex  hll^a  1  V  *  ^  spherically  blunted  cone  with  an 

TlL  1  ^  f  ^  °  «  =  10°  p  in  some  points  on 

com  f  ^  that  in  the  frozen  flow.  The 

abov^An  evaluating  acf.^al  flight  effects  at  altitudes 

above  60  km  and  at  velocities  of  6,  500  m/sec  or  greater. 
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Antonets,  A,  V, ,  A,  V,  Krasil'nikov,  and 
V,  I,  Lagutin.  Experimental  determination  of 
the  position  of  pressure  center  in  hypersonic 
gas  flow  around  a  blunted  cone  at  different 
angles  of  attack.  MZhiG,  no.  2,  1971,  142-143. 

Following  from  the  previous  paper,  the  relationship  between 
the  position  of  a  pressure  center  Cj  and  an  angle  of  attack  or  was  determined 
in  a  hypersonic  wind  tunnel  at  M«q  =  6  and  15  and  Reynolds  number  Re  BIO^ 
using  a  spherically  blunted  test  cone.  The  cone  had  a  sphere  radius  R, 
relative  conic  length  x/R  =  5,  and  apex  half-anglo  6=  11°.  The  highly  accurate 
method  of  free  aerodynamic  balance  was  applied  while  varying  the  distance 
Xq  of  the  rotational  axis  from  the  cone  nose.  The  Cj  value  for  a  balanced  a 
coincides  with  Xq.  Plotted  experimental  Cj  data  (Fig.  1)  agree  within  a  2% 


Fig.  1.  Pressure  center  versus  angle 
of  attack:  curves  1  and  2  -  calculated  for 
Mao  =  6  and  15,  respectively;  small  circles 
>  experimental  data  for  Mco=  6;  points  - 
experimental  data  for  Mgo  =  15. 


maximum  discrepancy  with. theoretical  data  calculated  by  the  method  of 
Antonets  cited  in  the  previous  paper  [MZhiG,  no.  2,  1970,  178],  but 
disagree  with  the  Cj  ( a )  data  calculated  by  the  Newton  formula  (dotted  line). 
Results  suggest  that  physical  phenomena  on  the  leeward  side  of  a  cone  at 
o  2  9  do  not  significantly  affect  aerodynamic  characteristics. 


Voronkin,  V.  G.  Nonequilibrium  viscous 

flow  of  a  multicomponent  gas  near  the  stagnation 

point  of  a  blunt  body.  MZhiG,  no.  2,  1971,  144-147. 

A  theoretical  study  is  made  of  the  effect  of  shock  layer  variable 
parameters  on  heat  transfer,  in  the  case  of  nonequilibrium-dissociated 
air  at  a  barrier  of  finite  catalytic  activity  near  the  stagnation  point  of  a 
spherically  blunted  body.  Heat  transfer  in  the  nonequilibrium  outer  flow  is 
compared  to  that  in  an  equilibrium-dissociated  gas  flow.  Numerical 
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calculation  of  the  nonequlllbrlum  flow  was  based  on  the  concept  of  a  compressed 
viscous  shock  layer.  Only  free-stream  Reynolds  numbers  R  above  1,500 
were  considered.  Under  these  assumptions,  the  calculation  was  reduced 
to  a  solution  of  the  boundary  layer  equations  with  boundary  conditions  corre¬ 
sponding  to  those  Immediately  behind  the  leading  shock  wave.  Temperature 
profll^Bs,  concentrations  of  oxygen  and  nitrogen  atoms  In  the  shock  layer, 
and  the  relative  heat  flux  were  calculated  at  a  60  km  height,  7.  4  km/sec 
free-stream  velocity,  700°K  wall  temperature  and  1-100  cm  radius  of 
curvature  R.  It  was  shown  that  In  the  nonequlllbrlum  shock  layer  at 
small  R  can  be  significantly  Increased  (by  a  factor  of  3)  In  relation  to  the 
theoretical  In  an  asymptotically  thin  nonequlllbrlum  boundary  layer  with 
equilibrium  parameters  at  Its  outer  boundary.  At  Increased  R  the  vortlclty 
effect  decreases  and  the  effects  of  nonequlllbrlum  chemical  reactions  become 
predominant. 


Lunev,  V.  V. ,  and  N.  Ye.  Khramov.  Flow 
past  a  blunt  body  near  the  stagnation  pclnt  in 
a  divergent  hypersonic  stream.  MZhlG,  no.  3. 
1970,  102-105. 


,  .  ^  Inhomogeneous  hypersonic  outer  stream  near  the  stagnation 

of  analyzed  at  a  small  ratio  k  of  densities  ahead 

of  and  behind  a  straight  shock  wave.  Thlckncsu  6  of  the  shock  layer  at  the 
forward  stagnation  point  Is  formulated  by 


6  1 


1+V!24- 


VSElO 

(v«l) 


(1) 


where  R*  Is  the  effective  curvature  radius  of  the  body  and  v  =  0  and  v  =  1 
refer  to  two-dimensional  and  axlsymmetric  flows.  The  formulas  In  (1) 
were  derived  frorn  equations  of  motion  assuming  gas  density  to  be  constant. 

It  follows  from  (1)  that  8,  and  hence  the  profile  of  the  shock  layer,  depend 
on  k  and  R^  only.  This  dependence  Is  Illustrated  by  computed  6(k)  plots 

strongly  underexpanded  conic  jet  around  a  sphere. 

The  l^nitlal  flow  parameters  are  tabulated.  The  data  all  He  on  t^e  same  curve 
which  coincides  with  the  theoretical  curve  calculated  for  v  =  1  and  R*/R 

x/R"r!t?  pressure  distribution  p/p^  along  a  barrier  versus 

x/R  ratio  diverge  slgnlflcmtly  for  different  jets,  but  the  same  plots  of  p/p' 
versus  x/R^  form  a  relatively  narrow  band  which  Is  close  to  the  plot  for  ° 
a  sphere  In  a  homogeneous  stream  and  the  Newtonian  plot.  The  effect  of 
stream  divergence  Is  thus  accounted  for  by  Introduction  of  R*.  The  effects 
of  nonuniform  density  and  velocity  of  the  oncoming  stream  are  taken  Into 
account  In  an  analogous  manner* 


/ 
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Arutyunyan,  G,  M,  Climb  of  a  ahock  wave  on 
a  wedge  moving  at  auperaonlc  velocity,  IN; 

Volny  V  neuprugikh  aredakh  (Wavei  in  Nonelaatic 
media),  Kiahinev,  AN  Mold  SSR,  1970,  7-12, 

The  flow  generated  at  the  upper  aurface  of  a  moving  wedge  by 
an  oncoming  plane  ahock  wave  waa  analyzed.  Conditiona  are  formulated  for 
the  flow  according  to  the  pattern  ahown  in  Fig,  1.  It  waa  aaaumed  that  intenaity 
of  the  oncoming  ahock  wave  CE  waa  equal  to  that  of  the  attached  compreaaion 
ahock  OA,  CE  ia  normal  to  OB,  and  the  interacting  ahocka  CD  and  DA  are 
weak  and  reflected  at  a  )^/2  angle  by  a  rigid  wall  with  reaulting  refraction, 
but  without  formation  of  a  third  ahock  wave.  It  ia  ahown  that,  under  theae 
aaaumptiona,  the  flow  paat  the  wedge  upper  aurface  OB  can  be  calculated  uaing 
linearization  theory,  if  a  point,  defined  by  the  Mach  number  M  of  the  oncoming 
flow  and  the  angle  of  the  wedge,  liea  within  certain  boundariea  in  the  M-6 
plane.  All  the  conditiona  were  verified  within  an  area  bound  by  the  curvea 

fl  (M,  8  )  =  0,  f2  (M,  8  )•  (M,  8  ).  and  the  abaciaaa.  If  the  point  defined 


Fig.  1.  Flow  diagram  on  the  upper 
aurface  of  a  wedge. 


by  the  wedge  parametera  M  and  8  liea  within  thia  area,  the  flow  will  be 
uniform.  Preaaure  within  the  area  ia  then  given  by 


P* 


2PM* 


(P 


(1) 


where  Po  i*  the  preaaure  generated  by  interacting  ahocka  and  Y  ia  the 
adiabatic  exponent. 
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Mamadaliyev,  N.  A. ,  and  Kh.  A.  Rakhmatufln. 
Streamline  fluid- aolid,  two-phase  flow  pait 
>  Blender  profile.  IN:  Volny  v  neuprugikh 
Bredakh  (Waves  in  nonelastic  mr  dia).  Kishinev, 
AN  Mold  SSR,  1970,  146-152. 


elenHor  n,  «i  of  a  fluid-solid  supersonic  flow  past  a 

oiramlter^  **  f®*-  “a®  in  calculations  of  jet 

parameters  from  supersonic  nozzles,  flight  mechanics  in  a  gas-dust  medium 

th^orlfne""  r  V  post-shock  reflection  of  particles  from  ' 

profile  surface.  For  a  specific  volume  <1,  the  effects  of  particles  on  the 

gas  now  and  friction  between  particles  in  the  region  of  a  three-velocity 

Da°rHcrr.  1*  Solution  of  linearized  equations  of  motion  (or 

particles  impinging  on  and  reHected  from  a  wedge  with  a  given  apex  angle 

a^d  ‘he  boundary  between  the  three- velocHy 

and  the  two- velocity  flows  is  a  straight  line  with  slope  =  2  R  The  ^ 

boundary  condition  at  y  =  0  is  also  a  straight  line,  and  formulaies  the  law 
of  mirror  reflection.  For  the  case  of  a  two-phase  medium  with  a  significant 

Tso  n  fort^rhe*?  *  ‘he  effect  of  the  reflected  stream, 

the  nrX  !  r  7“*  '**'°^**"‘  ‘he  form  of  an  equation  oi 

a^paration  boundary  with  a  constant 

inclination  .  In  contrast  to  the  first  analyzed  case,  the  profile  surface  is 
weakly  concave  in  relation  to  the  stream,  and  the  three-velocity  region  II 
narrows  with  the  increase  in  X|,  as  shown  in  Fig,  1, 


Fig.  1,  1  -  separation  line,  profile  surface; 
2  -  at  a  small  and,  3  -  at  a  significant 
particle  concentration. 
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Lunev,  V,  V,  Nonvlscoue  flow  noar  the 
critical  point  of  a  blunt  body  having  a  variable 
denaity  along  ita  axia.  MZhlG.  no.  A.  1Q71 
56-61, 

The  flow  in  a  thin  compreaaed  layer  around  a  amooth  blunt 
body  ia  deacribed,  in  the  framework  of  hyperaonic  theory.  A  aet  of 

equationa  ia  derived  uaing  an  x,  y  curvilinear  coordinate  ayatem 
with  the  atagnation  point  aa  origin.  A  localized  cloaed  aolution  near  the 
atagnatlon  point  ia  poaaible  when  the  ahape  of  the  ahock  wave  approachea 

w.'v.'.V'ir  ''  •"'1  behind  Ih.  .hock 

w.»e  lend,  toward,  zero.  A  .olullon  In  the  vicinity  of  the  .xl.  1.  .ought 
in  the  form  of  an  expanaion  in  xK  powera: 

~w.(y) r  — x  +  ... 

Ap-rp.(v)+... 

where  K  =  R-1  i,  the  body  curvature,  uV  and  vV  are  velocity  componenta 

aLa*  H***'.f!  ‘n*  now,  r  ia  the  diatance  from  the  aymmetry 

H  i  A  u."  U  ^  'Obititution  of  variablea,  conaervation  equationa 
«e  derived,  which  are  aimilar  to  adiabatic  flow  equationa  for  a  perfect  gaa 

Buaemann  formula  from  (1)  * 

Ind  cated  that  univeraality  of  preaaure  diatribution  can  be  extended  to  the 

♦I?.  !*  of  a  amooth  body.  The  Buaemann  formula  alao  led  to 

the  concluaion  that  the  aaymptotic  preaaure  diatribution  along  the  body  ia 

no.*.m  t,  A  Analyol.  of  boundary  condition,  revealed  two 

poealble  boundary  flow  pattern..  According  to  one  pattern,  .  dl.trlbullon  In 
the  vicinity  of  the  wall  aatl.flea  the  condition 

9  ^const  k  (2 ) 

*1^®  ■“^“non  in  aimilar  variablea  f ,  t,  including  t  6,  dependa  only  on 
the  body  valuea  of  a  diatribution  and  the  axial  preaaure  gradie^  (at  f  =  0) 

*■  the  aolution  in  the  caae  of  a  conatant  p  =  kj*  along  the 
Wckneaa  6  of  a  ahock  layere  The  aolution  of  (2)  for  an  axiaymmetri^l  flow 

If  the°ord  V  weakly  dependent  on  k„.  except  for  a  narrow  near- wall  region 
of  the  order  t  -  k  J  .  The  aolid  linoa  in  Fig.  1  are  plotted  for  the  condition 

r:o*j/"n%*r «>..ocUt.d’-.lr  .f.“n-.,umtZ  “ 

1  5  dene  d  ’  **’*  •’•ickne.e  6/kR  calculated  for  curve. 

1.5  depend.  .Ignlflcantly  on  ,  dlolrlbutlon.  The  .olullon  for  a  two-dlmen.lon.l 


Fig.  1,  Axial  velocity  profile!  versus 
initial  coordinate  y/6  ;  the  nunnbers  1-5 
at  the  curves  correspond  to  those  in  Fig.  2. 


Fig.  2,  Density  profiles  versus  initial 
coordinate  y/J  for  an  axisymmetric  flow. 


now  (v  =  0)  is  seen  to  be  always  dependent  on  k^;  but  the  dependence  on 
ko  J  0.05  -  0.2  is  weak  as  shown  in  Fig.  3,  where  the  plots  1-4  correspond 
o  o  =  0.2,  0.05,  0.01,  and  0.001.  A  certain  degree  of  universality  of  the 

*  J  /  *lft  )  function  is  consequently  also  to  be  expected  in  the  case  of  two 
-dimensional  How. 


Axial  velocity  profiles  versus 
similar  coordinates  f,  t,  for  a  two-dimensional 
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According  to  the  second  flow  pattern,  P  distribution  near  the 
wall  exhibits  a  singularity.  A  significant  variation  occurs  in  the  p  profile 
which  is  bound  to  affect  the  flow,  particularly  in  the  near>wall  layer.  Even 
i.n  this  case,  however,  the  f  ( tj  /»|  g  )  prc^jie  is  expected  to  remain  universal, 
because  of  the  weak  dependence  of  the  solution  on  the  parameter  k^,  i.e. , 
the  gas  p  near  the  wall. 

It  is  concluded  that  in  all  except  the  i.ear>wall  sublayer 
case,  the  solution  should  be  the  same  as  for  a  constant  p  equal  to  the  p 
value  behind  the  shock  wave.  Calculation  of  the  structure  of  an  inviscid 
shock  layer  can  thus  be  simplified,  given  the  universality  of  the  axial 
flow  profile  in  f,*!  variables. 


Neyland,  V.  Ya.  Upstream  propagation  of 
perturbations  from  interaction  of  a  hypersonic 
flow  with  a  boundary  layer.  MZhiG,  no.  4,  1970, 
40-49. 


Hypersonic  viscous  gas  flow  past  a  slender  body  is  analyzed  for 
weak,  moderate  and  strong  free  interactions.  In  the  case  of  a  weak  hyper¬ 
sonic  interaction,  perturbations  generated  in  the  vicinity  of  a  singular  point, 

(e.  g. ,  bottom  cut  or  angular  point),  are  seen  to  propagate  upstream  a  short 
distance  only.  Locally  nonviscous  flow  regions  may  also  form  with  high 
pressure  gradients.  Boundary  layer  separation  into  sub-regions  of  different 
sizes  disf-ppears  in  the  presence  of  a  moderate  or  strong  hypersonic  interaction. 
The  perturbations  in  this  case  propagate  up  to  the  leading  edge  of  the  body  in 
the  process  of  free  interaction  along  the  entire  boundary  layer.  Proof  was 
obtained  that  localized  flow  regions  of  very  high  pressure  gradients  are 

not  generated  at  the  interaction  parameters  X  i  1.  The  possibility  exists, 
however,  of  a  separation  zone  formation  on  the  order  of  magnitude  of  one 
body  length.  The  hypersonic  boundary  layer  at  X2  1  is  described  in  a  first 
approximation  by  a  set  of  equations  derived  from  Navier-Stokes  equations  * 
with  boundary  conditions.  Boundary  conditions  at  the  aft  end  of  a  body  must 
be  satisfied  in  addition  to  those  near  the  leading  edge.  A  complete  set  of 
boundary  conditions  was  established  in  the  form  of  a  nonsingular  solution. 

The  conclusion  was  drawn  that  the  flow  at  and  X  2  1  is  described 

by  the  same  boundary  layer  equations  along  the  entire  separation  zone, 
except  for  small  regions  such  as  those  near  the  angular  points. 
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Zakharchenko,  V.  F.  ,  M.  V.  Tsvetkova, 
and  Ye,  E,  Borovskiy.  Suy>ersonic  flow 
around  a  porous  cone.  PM,  no.  2,  1972, 

132-134, 

Supersonic  air  flow  around  a  circular  pointed  cone  with  a 
permeable  surface  was  investigated  in  connection  with  the  control  of 
incident  heat  flux  and  drag  on  streamlined  bodies.  Shadow  photographs 
were  obtained  of  the  unperturbed  flow  around  a  0. 197  rad  half  angle  cone 
at  0-0.24  values  of  the  injection  parameter  (p  V)i,  where  pi  and  Vi  are 
the  density  and  velocity  of  the  injected  gas.  The  Mach  number  M  oo  of  the 
incident  unperturbed  flow  was  2.9,  the  Reynolds  number  Re«,  /m  was 
2.6  X  10  ,  and  the  ratio  Tqoo  /Tqi  of  stagnation  temperatures  of  the  external 
and  injected  gas  was  1.  Photographs  were  used  to  determine  the  angle  0 
between  the  compression  shock  front  and  the  direction  of  the  unperturbed  ® 
flow.  It  was  assumed  that  the  free  stream  boundary  is  also  conic  at  uniform 
gas  injection  through  the  conic  surface  and  that  the  pressure  along  the 
boundary  is  uniform.  The  boundary  can  then  be  assimilated  to  an  impervious 
surface  of  an  imaginary  cone  as  well  as  its  position,  i.  e.  the  angle  A0 
between  the  boundary  and  the  cone  surface.  Pressure  p  along  the  boundary 
can  therefore  be  calculated  as  in  a  problem  of  supersonic  flow  around  a 
cone.  Calculations  indicate  that  both  0^  and  a0  increase  continuously  with 
increases  in  injection  (p  V)i,  the  rate  of  increase  depending  on  the  ( 
value.  The  pressure  coefficient  pg  on  the  boundary  was  calculated  using 
available  data  ot^supersonief  flow  around  cones.  Plots  of  the  relative  pressure 
coefficient  px  -  p/|f  along  the  cone  surface  at  different  (p  V)i  values  (Fig.  1) 
show  a  qualitative  pattern  of  variations  of  the  pressure  coefficient  in  the 
presence  of  air  injection  in  comparison  with  the  flow  (p)  around  an  impervious 
cone,  i. e,  in  the  absence  of  air  injection.  The  experimental  pj  -  p/p  value 
is  shown  (Fig,  2)  to  depend  significantly  on  the  coordinate  x/L,  and  (]rv)j. 


Fig.  1,  Relative  pressure  coefficient 
versus  injection  parameter. 


45 


Fig,  2,  Relative  pressure  coefficient 
versus  coordinate  x/L  at  fixed  points 
of  the  cone  generatrix. 


The  data  reveal  a  significant  dependence  of  incident  flow 
parameters  on  the  injection  pressure  and  make  it  possible  to  calculate  the 
drag  on  a  permeable  surface  cone.  Total  drag  calculations  should  include 
an  allowance  for  injection  drag  which  can  amount  to  10%  of  the  total  drag 
at  a  low  injection  pressure. 


LifshitSi  Yu,  B,  and  E,  G,  Shifrin,  On  the 
problem  of  transonic  flow  over  a  convex  corner, 

MZhiG,  no,  2,  1971,  67-69, 

Mathematical  analysis  of  two-dimensional  potential  flow  in  the 
vicinity  of  the  sonic  point  O  shows  that  a  similar  solution,  satisfying  the 
conditions  of  Vaglio- La urin's  solutionfj.  Fluid  Me ch, ,  i960,  v,  9,  no,  l], 
is  possible  for  equations  of  transonic  flow  over  a  convex  corner,  A  solution 
to  the  ^  (u,  v)  hodograph  plane  is  obtained  by  analogy  with  that  of  Tricomi 
equations  which  become  equal  to  zero  at  the  v  =  2/3  characteristic  of 

the  first  family  of  flow  lines.  A  solution  C  =  ( uj^/v^  to  the  0  equation 
was  found  in  the  third  quadrant  of  the  u,  v  plane.  It  is  concluded  that  the 
profile  in  the  subsonic  region  must  be  convex.  A  one- sheeted  plot  of  the 
solution  in  the  hodograph  plane  indicates  that  the  sonic  line  is  convex  in  relation 
to  the  region  of  subsonic  velocities.  The  surface  velocity  vector  at  a  transition 
point  from  the  subsonic  region  is  given  by  u  =  ,  where  K  is 

the  curvature  of  profile  in  the  O  point.  The  similar  solution  derived  satisfies 
a  boundary  condition  in  the  vicinity  of  the  salient  point,  as  well  as  the  condition 
of  finite  curvature  in  the  subsonic  region.  The  solution  also  describes  the 
principal  member  of  asymptotic  expansion  in  axisymmetric  and  rotational  flow. 
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I.  A.  Bogashchenko,  A.  V.  Gurevich,  R.  A. 

Salimov,  Yu,  I,  Eydelman.  Flow  of  a  rarefied 
plasma  around  a  body.  ZhETF,  v.  59,  no  T~ 

1970,  1540-1555. 

^  the  structure  of  the  disturbed  wake  behind  a  disc  around  which 
a  rarefied  plasma  is  flowing  is  studied  in  detail  by  the  Langmuir  probe 
^chnique.  The  plasma  was  produced  in  a  Q-machine  as  shown  in  Fig.  1. 
Properties  of  the  plasma  flow  are  discussed.  'A  characteristic  radial 
istribution  of  quasi-neutral  plasma  is  showh  to  exist,  with  maximal  density 
ai.'ound  the  system's  axis.  The  distribuflon  of  ions  in  the  plasma  is  given  by 


Fig.  1.  Plasma  generator. 

'■  1  -  ionizer;  2  -  anode;  3  -  disc  #1; 

'  ^  4  -  disc  #  2;  5  -  test  probe,' 

Maxwell's  function;  the  anpde  potential  is  chosen  so  that  the  electron  distri¬ 
bution  in  the  plasma  can  also  be  given  by  the  Maxwell  function.  The  density, 
velocity  and  temperature  in  the  plasma  flow,'  surfaee  potential  and  position  of 
the  disc  were  varied.  Test  measurements  and  their  comparison  v/ith  theory 
indicate  the  presence  of  ion  acceleration  in  a  selfconsistent  electric  field 
produced  in  the  plasma.  In  the  immediate  vicinity  behind  the  disc  the  effect 
of  the  electric  field  on  ion  motion  is  quite  pronounced,  causing  a  noticeable 
ion  concentration  near  the  axis.  The  plasma  wake  has  a  damped  oscillatory 

profile;  it  is  shown  that  the  decrease  of  ionizer  temperature  and  the  increase 
of  plasma  density  will  increase  the  damping,  and  for  Ts  1800OK  the  distribution 
of  plasma  behind  the  body  assumes  an  almost  monotoriic  character.  At  large 
distances  from  the  disc,  plasma  disturbance  is  shown  to  be  a  function  of  the 
radius  only.  An  important  role  in  oscillation  damping  is  played  by  the  defor¬ 
mation  of  the  ion  distribution  function.  The  effect  of  ionic  collisions  on  the 
ion  distribution  function  is  discussed  in  detail.  An  investigation  of  oscillation 
damping  shows  that  the  longitudinal  and  transverse  ion  temperatures  gradually 
level  out  as  a  result  of  collisions.  It  is  shown  that  the  change  of  sign  of  the 
body  potential  relative  to  that  of  the  plasma  does  not  appreciably  affect  the 
structure  of  the  disturbed  wake  behind  the  disc.  However,  a  change  of  sign 
in  the  potential  of  another  body  located  in  the  disturbed  zone  of  the  first  body 
strongly  affects  the  structure  of  the  zone  in  the  region  between  the  bodies. 
Theoretical  and  experimental  data  are  compared  for  the  axial  and  radial 
plasma  distribution  in  perturbed  regions,  and  the  results  of  the  measurements 
are  found  in  accordance  with  theory. 
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yil.  ABLATION  STUDIES 


Karasev,  A.  B. ,  and  T.  V.  Kondranin.  The 
effect  of  ablation  products  on  heat  transfer 
during  destruction  of  graphite  in  radiative 
air  plasma.  M^iG,  no.  1,  1971,  23-31. 


air  around  th.  l  ^  ^  hypersonic  flow  of  viscous,  heat- conductive,  radiating 
th^  t  t  of  a  graphite  sphere  is  analyzed  in 

Z  the  separated  shock  wave  and  the  body  surface.  The  flow 

s  described  by  a  set  of  equations  with  self- similar  variables  and  boundary 
radia  the  shock  wave  and  on  the  surface.  Heat  transfer  in  the 

annt  •  calculated  by  solving  the  equations  by  the  method  of  successive 

Chemfc^  adopted  for  diffusion  of  ablation  products. 

Chemmal  reactions  were  considered  to  be  frozen  in  the  boundary  layer  and 

equilibrated  at  the  boundary  layer- shock  wave  interface.  Gas  state  on  the 
ablation  surface  was  determined  from  the  chemical  equilibrium  condition, 
assuming  equality  of  vapor  pressure  with  stagnation  pressure.  Radiative 
and  convective  heat  transfer  were  calculated  in  the  presence  or  absence  of 
ablation.  Analysis  of  plotted  data  shows  that  the  spectral  range  of  radiant 

arhifh  Ifm  streamline  flow  around  an  impervious  surface 

at  high  atmospheric  reentry  velocities.  In  the  absence  of  ablation,  the 
t^otal  radiant  heat  flux  changes  insignificantly  by  transfer  through  a  cold 
canabf''^f  T'’*  presence  at  the  boundary  layer  of  ablation  products 

wHMncrl  radiation  means  that  heat  flow  to  the  surface 

11  increase  if  the  vapor  product  is  opaque  in  the  visible  range,  but  will 

^  Actually,  the  shielding 

n  Vi  s  ®  cancelled  by  the  radiative  effect  at 

nuMt.l,  d  by  these  authors  on  graphite  ablation  was 

publishtd  in  tiffects  of  Strong  Explosions,  no.  1,  1971,  p.  41, 


Karasev,  A.  B. ,  and  A.  N.  Lyakh.  Study 
of  radiative  and  convective  heat  transfer 
by  an  emissive  mixture  of  carbon  dioxide  and 
nitrogen  flowing  past  the  stagnation  point. 
MZhiG,  no.  2,  1971,  39-47. 


Hypersonic  flow  characteristics  of  a  viscous  thermallv 
conductive  mixture  (90%  COj,  10%  N2)  around  a  spherically  blunted  body 

1  •'"u,  radiation  energy  transfer  and  reabsorption. 

Self-similar  variables  and  boundary  condition  equations  describe  the  gas 
low  in  a  compressed  thin  shock  layer,  assuming  local  thermodynamic 
equilibrium  through  the  entire  shock  layer.  A  solution  was  obtained  by 
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iteration  and  successive  computation  for  flights  of  the  sphere;-  at  9,000  - 
12,000°K  and  l-lO  atm  pressure  behind  the  shock  wave,  -The  spheres 
were  of  1  m  radius  and  had  a  3000°K  wall  temperature,  Analysi,s  of  resulting 
data  revealed  an  insignificant  radiation  effect  on  convective  heat  transfer. 

The  relative  contribution  of  different  spectral  ranges  to  the  radiative  heat 
flux  on  the  wall  was  found  to  be  approximately  43,  40,  20,  and  15%  for  the 
0,  128  -  0,33^  (CO  band),  0,33  -  0,  66/-(  (CN  bands),  0,80  -  l,15/<and 
0,45  -  0,80^  ranges,  respectively.  The  contribution  of  the  UV  spectral 
range  decreases  while  that  of  the  visible  range  increases  with  an  increase 
in  pressure.  Radiation  heating  is  five  times  greater  than  convection  heating 
at  lO  atm  and  12,000OK,  Convective  and  radiative  heat  fluxes  in  the  CO2-N2 
mixture  exceed  the  corresponding  values  in  air  by  20  and  40%, 


Apshteyn,  E,  Z, ,  L,  G,  Tefimova,  and 
G,  A,  Tirsky.  Intensive  destruction  of  a 
vitreous  body  by  Cshock3  radiation.  MZhiG, 
no,  2,  1971,  131-134, 

Destruction  of  a  vitreous  body  by  radiation  from  a  shock  layer 
in  hypersonic  gas  flow  is  analyzed  assuming  an  optically  thin  boundary  layer 
at  high  injection  rates  and  sufficiently  intense  vaporization.  Presence  of  a 
liquid  film  near  the  stagnation  point  on  the  body  surface  is  presumed,  to 
simplify  calculations  of  destruction  parameters.  Application  of  an  approximation 
of  radiative  heat  conduction  in  the  body  to  the  liquid  film  led  to  a  set  of 
transcendental  equations  which  describe  ablation,  friction,  and  surface 
temperature  as  functions  of  shock  layer  parameters.  Iterative  solution  of  the 
equations  is  illustrated  by  diagrams  of  the  destruction  rate  and  gasification 
r  of  a  heat  insulating  material  with  thermophysical  characteristics  of  fused 
quartz  under  assumed  flight  conditions.  Plots  of  Fversus  body  radius 
Rg  show  that  F decreases  with  a  decrease  in  due  to  liquid  phase  ablation. 

At  great  Rg  values,  coating  destruction  proceeds  by  vaporization  only.  The 
contribution  of  liquid  phase  ablation  to  the  destruction  rate  becomes  significant 
at  small  Rg,  because  of  the  effect  of  the  radiative  heat  conductance  parameter 
^  R  Since  it  increases  greatly  with  increase  in  atmospheric  reentry 

velocity,  radiative  heat  flux  is  thus  a  decisive  factor  in  calculation  of 
destruction.  A  previous  study  on  this  problem  was  reported  by  Apshteyn 
in  Effects  of  Strong  Explosions,  no.  1,  1971,  p.  22. 
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Vlll.  HIGH-PRESSURE  EFFECTS 

I 


Vereshchagin,  L.  F. ,  A,  A,  Semerchan, 

N,  N,  Kuzin,  and  Yu.  A,  Sadkov.  Performance 
of  a  three- stage  high  pressure  unit  with  an 
operational  volume  of  approximately  100  cm^. 

DAN  SSSR,  V.  202,  no.  1,  1972,  74-75. 

The  performance  characteristics  are  discussed  of  a  high 
pressure,  three  stage  unit  designed  to  produce  large  blocks  of  polycrystalline 
diamond  materials  and  composites  based  on  synthetic  polycrystalline  diamonds. 
These  materials  are  suggested  for  the  tooling  of  the  main  components  of 
another  high  pressure  unit  (with  an  operational  volume  of  approximately 
0,25  cm  ),  The  use  of  these  materials,  with  superior  physical  and  mechanical 
properties,  will  lead  to  a  substantial  increase  in  the  pressure  range  of  the 
0.  25  cm^  unit. 

Pressures  of  up  to  lOO  kbar  can  be  realized  with  the  100  cm^ 
unit.  Calibration  of  the  unit  pressure  scale  was  done  using  the  reference 
points  of  Bij.!!,  Bajj_jjj,  Biju.iv,  and  Sn;  and  also  by  the  method  of 
electrical  resistance  using  a  10,000  ton  press.  Concavity  of  the  pressure 
calibration  curve  is  caused  by  a  relatively  small  ratio  (2:1)  of  F/Fj,  where 
F  is  the  total  force  and  Fj  is  the  effective  first  stage  force,  (F/Fj  ratio 
is  as  much  as  5:1  for  the  0.25  cm^  high  pressure  unit. )  The  increase  of  the 
pressure  range  and  the  concavity  of  the  pressure  calibration  curve  are 
dependent  on  an  increase  of  the  F/Fj  ratio,  which  is  not  attainable  for  the 
loo  cm  unit  without  a  substantial  increase  in  the  second  and  third  stage 
cross-sections.  Since  larger  parts  made  from  the  presently  used  cermet' 
hard  alloy  (VK-6)  would  likely  suffer  a  loss  in  quality  due  to  sintering,  the 
effective  development  of  multistage  units  is  predicated  on  fabricating  at 
least  one  of  the  stages  using  high  strength  steel. 

The  authors  note  that  they  have  successfully  produced  large 
blocks  of  poly  crystalline  diamond  materials  of  the  ballas  and  carbonado 
types  using  the  100  cm^  high  pressure  unit. 
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Vereshchagin,  L.  F, ,  E.  Ya,  Atabayeva,  and  N.  A.  Bendeliani. 
Phase  diagram  of  Bi?Te'^  at  high  pressures  and  tenriperatures. 
FTT,  no.  8,  1971,  2452-2454. 


The  phase  diagram  of  semiconducting  Bi2Te3  in  the  40-85  kbar  and 
20-700°C  ranges  was  plotted  on  the  basis  of  the  experimental  isobars  and 
isotherms  of  electrical  resistance  R.  The  R(P)  and  R(T)  functions  were 
measured  mostly  on  single-crystal  p  -  or  n-type  Bi2Te3  wafers  in  a  high- 
pressure  chamber  designed  at  the  Academy's  Institute  of  High-Pressure  Physics. 
Lithographic  stone  or  AgCl  were  used  as  the  pressure-transmitting  medium. 

The  boundaries  of  three  high-pressure  phases  are  shown  on  the  P-T  diagram 
of  Fig.  1. 


Fig.  1.  Phase  diagram  of  Bi2Te3. 

1  -  4  -  electrical  resistance  method, 

5  -  DTA  method,  6  -  data  of  D.L.  Ball 
(Inorg.  chem. ,  1.  4,  1962) 

The  phase  transition  curves  are  described  by  the  equations 

/*!_„  =  68.5  —  0.04937',  =  82.8  —  0.0703 T*.  (1) 

=  84.0  —  0.0726  T ,  =  56.1  —  0.0094  7*, 

and  the  P-T  coordinates  of  triple  points  are  tabulated  (see  Table  I  below). 


Equilibrium 

^•kbar 

J  r,  •€ 

1 

Llq.+I+ii 

40.4 

570 

Llq.+ii+ni 

43.0 

560 

Ll^.+ill+iv 

50,8 

553 

II+III+IV 

52.0 

433 

Table  1. 


V 


52 


The  I  -  II  and  II  -  III  traniitions  are  reveraible;  the  II  -  IV  traniition  ii 
not  reversible  below  ISQOC.  On  cooling,  the  metallic  phase  IV  yields  a 
metastable  phase  at  20°C  and  25  kbar  pressure,  intermediate  between  the 
phase  IV  and  the  original  phase.  The  phase  transitibns  detected  are  of  the 
second  order. 


Kalitkin.  N.  N.  and  L.  V.  Kuz'mina.  Cold  compression 
curves  at  high  pressures.  FTT.  no.  8.1971.  2314-2318. 


^  f  E  of  elastic  compression  of  an  atom,  chemical 

potential,  electron  density  p.  and  intra-atomic  potential  are  calculated  at  0° 
for  each  element  of  the  Periodic  Table.  Calculations  based  on  a  quantum 

r°*  simpler  than  the  earlier  calculations 

from  ti^  electr^lc  spectra,  and  also  more  accurate  than  calculations  based 
on  the  Thomas-Fermi  model.  The  pressure  P.  as  the  other  thermodynamic 
quantities  cited,  was  expressed  as  a  function  of  p.  e.g. 


radius  of  the  Wigner-Seltz  cell  In  the  QSM.  The  tabulated  P 
and  E  data  for  13  elements  and  for  atomic  volumes  v  in  the  .100  +  0  to  .422  +  2 
range  wore  then  calculated  on  a  BESM-6  computer  >^ith  a  relative  error  in 
the  -3  to  +12%  range.  Analysis  of  the  calculated  numerical  data  with  allowance 
for  p  (r)  5:  Z(v'  at  ultra-high  pressures  led  to  an  interpolation  formula  P(Z.R) 
which,  together  with  a  formula  analogous  to  (3)  for  a  locallzable  P.  were 

""  ""<1  compound. 

can  be  plotted  with  the  help  of  the  cited  formulas  for  P  by  solving  the  set  of 
equations 


P(Z^,  W,)=aP(Z„  wj,  Ar>  1; 


(2) 


which  describe  the  equality  of  P  in  the  ceU  of  various  components  and  the 
volume  additivity,  respectively. 


Popov.  V.  N.  and  V.  I.  Chernyshov.  Equation 
of  state  of  a  helium-nitrogen  mixture.  IAN 
Energ.  no.  3.  1971.  162-165. 


A  procedure  for  statistical  mechanical  treatment  of  exDerimonfai 

Nz,  and  He-N2  mixtures  In  the  273  -  7730k  range  of  temperatures  and  at 
pressures  to  400  bar.  The  virlal  coefficients  B.  C.  D.  E  are  given  in  polv 
nomlal  form  as  functions  of  T.  The  coefficients  of  the  polynLlals  are 
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tabulated.  The  virial  equations  of  state  for  He.  N2.  and  He>N2  mixtures 
were  developed  from  the  described  treatment  of  the  experimental  data.  The 
equation  of  state  for  a  He>N2  mixture  where  the  index  m  refers  to  the  mixture, 

»m-  1  +  "rrf*  h  I-  +  r^n*.  ( 1 ) 

describes  the  experimental  results  within  +  0.20%.  Compressibility  factors 
and  density  p  values  calculated  from  (1)  were  in  catlsfactory  agreement 
with  earlier  experimental  values  over  the  experimental  temperature  and 
pressure  ranges.  This  fact  indicates  that  the  caloric  data  calculated  from 
(1)  are  sufficiently  accurate  and  reliable  for  design  of  technological  processes 
and  solution  of  structural  design  problems. 


Polyakov,  L.  M.  Effect  of  uniform  pressure 
on  the  plasticity  of  rock  salt  crystals.  UFZh, 

V.  17,  no.  2,  1972,  210-217. 

Results  are  described  of  tests  on  the  effects  of  uniform  pressure 
on  plastic  deformation  of  rock  salt  crystals.  The  tests  correlated  pressure 
with  the  appearance  of  defects  affecting  dispersion  and  deflection  of  light, 
and  also  with  threshold  and  type  of  destruction.  An  experiment  was  conducted 
with  natural  optically  homogeneous  rock  salt  crystals  c*  &  x  6  x  12  mm. 
Measurements  of  the  intensity  of  light  dispersion,  i  corresponding  to  the 
plastic  deformation  at  given  values  of  uniform  pressure  showed  a  peculiar 
process  of  destruction  of  NaCl  crystals.  In  the  beginning  stage  deformation 
takes  place  without  a  significant  change  in  light  dispersion;  then  at  some 
deformation  <ip.  intensity  starts  increasing  and  at  some  further  deformation 
^2P'  )p  curve  reaches  a  saturation  region,  which  finally  concludes  in 

the  ultimate  destruction  of  the  crystal  specimen  at  deformation  C3p.  The 
relationship  between  intensity  of  light  dispersion  and  extent  of  deformation 
is  shown  in  Fig.  1  with  increase  in  uniform  pressure,  the  intensity  of  light 
dispersion  decreases,  and  its  distribution  in  the  crystal  becomes  homogeneous. 


Fig.  1.  Relationship  between  intensity 
of  light  dispersion,  i  and  deformation  stage  C: 

1-0;  2  -  500  kg/cm2;  3  -  1000  kg/cm2; 

4  -  2000  kg/cm^;  5  -  3000  kg/cm^;  6  -  5  000  kg/cm^. 


def  H?  «  ?  coincides  with  the  formation  of  crystal 

Haht  c-  250  ^  ‘^Onlflcantly  lest  than  the  wavelength  of 

®  uniaxial  compression,  dispersion 
^ looo^io  /  •^‘‘e^rmatlon  level  £j  -0.65%;  if  deformed  under  pressure 
0  3000  kg,  cm  ,  then  deformation  Increases  to  ftp  »  8.9%.  Mathematical 

.T:*,,'.  r,‘  ■nnu.„c*.%f  preLr.  on  pr«“cuy 

«nd  the  value  o(  deformation  corresponding  to  light  dispersion.  The 

pl^lenwdT  ‘n  tletall  and  microphotographs  are 


Fateyeva,  N,  5,,  and  L.  F.  Vereshchagin.  Meltinc 
curve  of  molybdenum  up  to  90  kbar.  ZhETF  P,  v.  14 
no.  4,  1971,  233-235. 


A  teat  la  brieny  deacrlbed  In  which  the  melting  curve  of  pure  Mo  unde 
high  preaaure  waa  measured  by  an  optical  method  described  earlier  by  the 
authors.  The  method  is  based  on  simultaneous  determinatior.  of  the  radiation 
intensity  ratios  I1/I2  and  12/13  of  two  pairs  of  narrow  spectral  regions,  and 
their  subsequent  comparison  with  Planck's  law.  The  experimental  T(P)  plots 
of  Fig,  1  can  be  presented  by  the  linear  equation 


^  ■  2BI3  ♦  0,|«*]0-Sp 


(1) 


where  T  is  -nelting  point  in  OR  and  P  is  the  pressure  in  bars.  The  probable 
error  is  +%  for  both  T  and  P  measurements.  A  similar  experiment  was 
reported  by  Fatayeva  et  al  on  melting  characteristics  of  graphite  under  high 
pressure  [Explosion  Effects  Report  No.  2,  p.  58].  ^ 


tto  ^kbar 


Fig.  1.  Melting  curve  of  Mo  up  to  90  kbar, 
calculated  from  all  experimental  points. 

+  -  temperature  data  of  Ii/r2  measurements, 
X  -  temperature  data  of  I2/I3  measurements. 
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IX.  ELECTROMAGNETIC  FIELD-PLASMA  INTIRACTION 


Liberman,  M,  A. ,  and  A.  T,  Rakhimov. 

Penctrntion  of  c- m  waves  into  a  plasma  with 
allowance  for  nonlinearity.  ZhETF,  v.  61, 
no.  3,  1971,  1047-1056. 

The  authors  examine  the  structure  of  an  alternatihg  electro¬ 
magnetic  field  interacting  with  a  weakly  ionieed  plasma,  tox  the  case  of 
thermal  nonequilibrium  in  the  plasma.  It  is  noted  that  nonlinear  effects 
begin  to  appear  in  such  a  case  at  relatively  weak  fields  in  comparison  to 
the  characteristic  plasma  field.  The  discussion  is  limited  also  to  the  case 
where  the  incident  e-m  field  frequency  is  well  above  the  electron-atom 
collision  frequency,  but  below  plasma  frequency.  Using  this  model  the 
authors  show  that  the  effect  of  the  e-m  field  on  the  local  ion-recombination 
balance  in  the  plasma  causes  an  appreciable  change  in  the  penetration 
depth  of  the  field  into  the  plasma,  according  to  the  relation  c/cop  where  (i> 
is  electron  frequency  in  the  interaction  region.  Expressions  are  derived^ 
illustrating  the  nonlinear  decrement  in  e-m  intensity  with  penetration,  and 
some  sample  calculations  are  given  for  the  effect,  using  typical  plasma 
and  beam  parameters. 


Colant,  V.  Ye.  Wave  penetration  in  plasma 
at  frequencies  npar  the  lower  hybrid.  ZhTF, 
no.  12,  1971,  2492‘2503. 

Theoretical  considerations  are  presented  for  optimizing  the 
introduction  of  e-m  radiation  into  a  magnetized  plasma.  It  has  been 
shown  that  injection  near  the  lower  hybrid  frequencies  may  be  advantageous 
since  plasma  opacity  to  the  incident  wave  will  be  minimal  here,  to  a  wave 
with  correctly  applied  longitudinal  delay.  The  author  therefore  investigates 
the  transformation  region  for  frequencies  near  the  lower  hybrid,  an 
derives  expressions  for  optimum  energy  transfer  under  these  conditions. 

It  is  assumed  that  the  applied  frequency  lies  well  below  electron  cyclotron 
resonance  and  well  above  ion  cyclotron  resonance.  For  the  case  of  a 
uniform  applied  magnetic  field,  expressions  for  optimal  delay  structure 
parameters  are  obtained,  assuming  the  idealized  delay  configuration 
of  Fig.  1.  In  some  cases  in  a  uniform  field  the  requirements  on  delay 
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Figure  1,  Delay  geometry. 

1  -  delay  element;  2  -  vacuum; 

3  >  plasma. 

parameters  may  become  excessive;  in  that  event  a  nonuniform  external 
field  may  be  applied,  resulting  in  two  transfori.iation  regions  and 
simplifying  the  delay  line  operation. 
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X.  IONOSPHERIC  EFFECTS 


Bykov,  B,  P, ,  and  S,  I,  Kozlov.  Ion  and 
electron  kinetica  in  a  disturbed  ionosphere 
at  altitudes  of  100^200  km.  Geomagnetizm 
i  aeronomiya,  no.  2,  1972,  340-342. 

Results  are  described  of  an  investigation  on  kinetics  of  ions 
and  electrons  at  altitudes  lOO  <  h  <  200  km  in  a  disturbr  1  atmosphere.  The 
disturbing  action  is  assumed  to  significantly  increase  the  degree  of  ionization 
O  of  air,  so  that  ©2^  and  O’*"  are  formed  approximately  in  direct 

relation  to  the  N2,  O2  and  O  contents  at  these  altitudes.  Investigation  is 
limited  to  values  a  1,  so  it  is  considered  that  the  atmosphere  is  not  warmed 
up  due  to  a  sharp  rise  of  a^i.  e.  electrons,  neutrons  and  ions  of  gas  have  the 
high  level  of  temperature  distribution  existing  in  normal  ionosphere. 


Form. 

no. 

Reaction 

Velocity  const,  cm^  per  sec 

1 

0+  +  Nj  KO*  +  N 

oi  =  2.10-» 

2 

0+  4.  Oj  -  +  0 

oj  =  -|.10-» 

3 

NV  +  O-O^  +  N, 

aj<  1.10-» 

4 

Nr  +0-»  N0++  N 

a»  «2,5.10-"> 

5 

n,*  +  Oj  -  0,’  4  Nj 

as  =  1 .  lO-w 

6 

Oj^  4  Nj  -  N0+  4  NO 

f:.< 

7 

0i+4«-»040 

a*  (0 .+)  =  2 . 2 . 1 0-f  (7’c/300)-> 

8 

N,-4«-N4N 

a’  (N,^)  =.  3- 10-’  (?V300)*'''* 

9 

NO-^  4  «  N  4  0 

o*  (NO^)  =  4.4. 10-’  (T,IZ00)-i,  200'- K 
<  7’,  <  500’  K 

a‘(NO+)  =  0.9. 10-»  (7’e/lOOO)-' T,  >  SOO”  K 

Table  1. 


The  general  scheme  of  elementary  processes  used  in  this  work 
and  their  velocity  constants  are  given  in  Table  1.  On  the  basis  of  these 
reactions  kinetic  differential  equations  are  obtained  which  were  solved  by 
computer.  Calculations  are  done  for  a  wide  range  of  altitudes  and  values  of 
^eo*  night  conditions.  Results  of  calculations  for  daylight 

conditions  at  03  =  lO”^^  and  =  10“^^  cm^/sec  are  plotted.  The  main 
characteristic  features  of  ion  and  electron  kinetics  in  a  disturbed  ionosphere 
are  seen  to  be  as  follows:  1)  The  role  of  N2^  ions  is  significant  only  at  the 
initial  moment  after  disturbance,  not  exceeding  in  most  cases  '-I  sec; 
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2)  At  lov/er  altitudes,  as  in  a  normal  atmosphere,  recombination  takes 
place  mainly  according*  to  a  square  law.  However,  with  an  increase  of 
Ngo  after  the  disappearance  of  over  a  certain  time,  velocity  variation 
may  follow  a  linear  law.  The  duration  of  this  time  significantly  increases 
at  higher  altitudes;  3)  The  time  behavior  of  fHO**"]!  d-iffers  little  from 
predictions.  Even  at  upper  altitude  regions,  the  role  of  this  ion  is  small, 
except  when  10^  cm“3;  4)  To  a  first  approximation,  at  all  altitudes 

the  main  molecular  ion  is  oj  :  5)  The  principal  difference  in  behavior  of 
ions  and  electrons  in  night  conditions  from  day  consists  in  the  fact  that 
the  dissociative  recombination  level  increases,  so  that  the  role  of  N2'*’  and 
NO**"  ions  substantially  decreases;  and  6)  The  decay  time  of  electron  concen¬ 
tration  from  N0O  to  normal  values  of  an  undisturbed  ionosphere  has  only  a 
slight  dependence  on  tiie  value  of 

These  conclusions  also  hold  up  under  varying  values  of  the 
unknown  velocity  constants  03  and  Ot^.  All  calculations  were  done  assuming 
Oj  -  to  be  independent  of  temperature. 


Alimov,  V.  A.  Effect  of  strong  radiowaves  on  the 
ionosphere.  Geomagnetizm  i  aeronomiya,  no.  I, 

1972,  346-348. 

The  author  proposes  possible  mechanisms  for  increasing  (or  decreasing) 
electron  density  in  a  localized  region  of  the  ionosphere  by  means  of  a  powerful 
r-f  beam,  and  shows  the  considerable  effect  that  such  a  beam  may  have  on  the 
nonuniform  composition  of  the  affected  region.  Jhe  analysis  gives  examples  of 
the  effects  obtainable  in  terms  of  the  D,  E  and  F  layers. 

As  is  known,  raising  the  electron  temperature  in  a  plasma  may  result 
in  a  drop  in  electron  density  N,  owing  to  ejection  of  plasma  from  the  heated  region, 
or  alternatively  N  may  increase  owing  to  an  upset  of  the  ion  recombination  balance. 
In  the  first  case  N  is  established  over  a  time  tD~ro/Di,  where  r^  is  the  distanc 
over  which  the  excitation  wave  amplitude  changes  significantly  [sicj ,  ^Jid  Di  is 
the  ion  diffusion  coefficient.  In  the  second  case,  N  is  determined  by  the  li  etime 
of  a  free  electron,  te  ~  N„)-^  where  a'  is  the  effective  coefficient  of  electron 

recombination  and  Nq  =  original  unexcited  density. 

Therefore,  for  to<  t^  the  plasma  ejection  process  is  dominant,  e.  g. 
for  the  upper  portion  of  the  F- layer,  wh^re 

N  A’i(l  -!-  (/?(r)  /A', 


(1) 
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in  which  E  =  amplitude  of  the  excitation  wave  and  Ep  =  plasma  field,  as  defined 
by  Gurevich[l],  Conversely,  when  te  <  tj)  the  recombination  process  dominates, 
as  in  the  lower  ionosphere  where 

=  f  (2) 

Using  standard  parameters  for  the  F-layer  and  assuming  r^  <  300  km, 
Alimov  shows  that  Eq.  (1)  applies;  for  the  E-layer  at  ro>  0.  1--1  km,  Eq.  (2) 
would  apply.  In  the  latter  case  it  must  be  emphasized  that  (2)  is  valid  only  for 
a  concrete  relation  of  a'  to  electron  temperature  Tg.  This  applies  in  the  E-layer, 
but  at  present  it  is  not  possible  to  assert  a  definite  a'  (Tg)  for  the  D-layer., 

A  brief  consideration  of  the  effect  of  ionospheric  drift  on  excitation 
lifetime  shows  that  drift  effect  is  negligible  on  the  steady- state  component  of  the 
excited  region,  providing  the  excitation  rate  is  sufficiently  larger  than  normal 
decay.  This  is  the  case,  for  example,  in  the  E>  and  F-layers  for  rg  >  7  km. 

Besides  changing  the  gross  level  of  N  in  the  foregoing  manner,  the 
r-f  excitation  also  affects  the  fluctuating  component.  Thus  in  the  E-layer  the 
mean  square  fluctuation  in  electron  density  will  increase  by 

rrN:;)Mi  -t- 


and  in  the  F-layer  will  decrease  by 


(A.Y)=  (A:Yo)v’ (1  H  (4) 

The  turbidity  of  the  excited  region  is  then  shown  to  be  a  function  of  r-f  field, 
where  turbidity  factor  a  is  given  by 

u  =  ±=.-.l_cxp(-j=).  (5) 

in  which  is  incident  energy,  is  scattered  energy,  and  is  mean  square 
phase  fluctuation  in  the  layer,  _Using  a.n  earlier  notationf2jfor  transmission  of 
r-f  through  the  layer,  namely  s2  =  ( A  N)2/(  A  N)2,  the  author  obtains 

graphical  solutions  from  (3),  (4)  and  (5),  showing  the  effect  of  applied  field  on 
turbidity.  These  are  given  in  Fig,  1  for  the  E-  and  F-layers  respectively, 
and  show  the  considerable  effect  of  incident  r-f  field.  This  ar/pears  particularly 
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f  ^  rtf  y  ;  j  i 

Fig.  1.  Turbidity  factor  variation  with 
applied  field  for  E-layer  (a)  and  F-layer  (b). 

significant  for  the  F-layer,  since  there  the  turbidity  effect  is  sharply  reduced  at 
sufficient  transmitted  power,  tending  to  nullify  the  disturbance  effect. 


References 


1.  Gurevich,  A,  V,  Geomagnetizm  i  aeronomiya,  v.  5,  1965,  70; 
V.  7,  1967,  291. 

2.  Alimov,  V.  A.,  et  al.  ,  ibid.,  v.  11,  1971,  790. 


XI.  RADIATION  DAMAGE 


Brudnyy,  V.  N,  ,  V.  P.  Voronkov,  M.  A.  Krivov, 
and  S.  V,  Malyanov,  Effect  of  electron  radiation 
on  gallium  arsenide  photodiodes.  IVUZ  Fiz,  no.  1, 
1972,  106-107. 


Results  are  reported  of  an  investigation  on  characteristics  of  gallium 
arsenide  photodiodes  after  radiation  by  electrons  with  1.5  MeV  at  SOOoK.  Experi¬ 
ments  were  done  with  diodes  having  a  base  thickness  of  120,  130  and  520  microns 
at  radiating  flux  densities  of  0;  1.  1  x  lOlScm'^,  2.  2  x  lO^^cm-Z  and  3.  3  x  10^5j,^-2, 
In  all  cases  the  volt-ampere,  spectral,  and  luminescence  characteristics  were 
graphed  for  conditions  before  and  after  irradiation  (four  characteristic  curves 
are  given).  Comparison  of  the  curves  shows  a  significant  change  in  the  characteristic, 
of  photodiodes  after  radiation,  particularly  in  the  forward  v~a  characteristic. 

The  authors  recommend  taking  these  effects  into  account  during  operation  of  such 
devices  in  a  radiation  field. 


Vorob  yev,  A.  A.,  0.  B.  Yevdokimov,  and  N.  P.  TubaJov, 
Tomsk  Poly technical  Institute  Imenl  S.  M.  Kirov. 

Effect  of  Overcharge  of  a  Dielectric  Charge  by  an 
Electron  Beam.  Leningrad,  FTT,  no.  12,  Dec.  1971. 
3691-3692.  ’ 


It  Is  shovm  for  the  first  time  that  a  change  In  the  sign  of  the 
effective  space  charge  can  occur  In  plexiglass  In  the  course  of  time.  It  Is 
noted  that  It  was  known  that  a  negative  apace  charge  of  thermally  Insulated 
electrons  Is  produced  upon  the  Irradiation  of  hlgh-ohmlc  dielectrics  by  fast 
electrons,  and  this  appears  In  the  development  of  discharge  diagrams  and  Is 
applied  in  studying  the  electrical  properties  of  dielectrics.  Plexiglass 
discs  4.5  mm  thick  and  coated  with  5-micron  aluminum  foil  were  irradiated  by 
0.8-1.2  Mev  electrons  at  a  current  density  of  0.5  pamp/cm2  for  15-20  sec,  and 
then  the  space  charge  was  tested  using  the  effect  of  an  electric  field  on  the 
passage  of  beta  particles  through  a  dielectric.  The  charge  samples  were 
periodically  irradiated  by  beta  particles  from  a  strontium-iridium  source, 
and  the  n^ber  of  passing  beta  particles  was  measured.  The  figure  is  given 
showing  the  relative  change  in  the  number  of  passing  beta  particles  due  to 
the  electric  field  as  a  function  of  time  after  irradiation.  The  figure  shows 
that  the  change  in  sign  of  the  effective  space  charge  is  observed  if  the 
maximum  in  the  distribution  of  the  stopped  electrons  shifts  to  one  of  the 
sur  aces  of  the  plate.  In  this  case  a  positive  space  charge  arises  in  the 
region  of  the  dielectric  not  occupied  by  the  negative  space  charge  due  to 
njection  from  the  surface.  Overcharge  is  attributed  to  leakage  of  the  nega- 
_^ve  space  charge  more  rapidly  than  leakage  of  the  positive  space  charge. 

The  observed  overcharge  effect  is  related  to  the  electric  state  of  the  plexi- 
glass,  just  as  in  polystyrene,  which,  as  distinct  from  plexiglass,  does  not 
yield  an  electric  state  and  overcharge  is  not  observed^ 
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XII.  LASER  EFFECTS  ON  MATERIALS 


Avotin,  S.  S.,  E.  P.  Krivchikova,  I.  I.  Papirov, 

P.  I.  Stoyev,  and  V.  I.  Tereshin.  Change  in 
electrical  resistance  of  berylliurm  from  laser 
radiation.  ZhETF,  v.  62,  no.  1,  1972,  288-293. 

Transient  and  long-term  effects  of  laser  irradiation  on  resistivity 
of  beryllium  are  described.  Tests  were  run  on  0.2  mm  thick  Be  foil 
strips,  exposed  to  pulsed  radiation  from  a  GOR-lOOM  laser  (not  identified) 
generating  1  millisecond  pulses  at  40  j,  focused  to  a  2.  5  mm  dia.  spot 
on  the  foil  surface.  Resistivity  p  was  measured  during  and  after  exposure 
with  a  constant  1  a  current  passed  through  the  specimen.  Variation  of  p 
during  exposure  is  illustrated  in  Fig.  1,  with  calculated  values  of  local 


Fig.  1.  Resistivity  variation  in  laser-irradiated 
beryllium. 

temperature  rise  ATgff  also  included.  A  detailed  examination  shows  three 
distinguishable  stages  in  A  p  (t)  up  to  the  TOO^ttsec  point,  as  indicated  by 
the  curve  form.  In  the  initial  stage  ( <  IOOjls)  a  local  heating  rate  of  about 
lO  deg/sec  was  calculated  to  exist  in  the  impact  region.  It  is  noteworthy 
that  a  full  anneal  back  to  pre-exposure  values  of  p  occurs  when  the  test  is 
done  :  -  a  SOQoK  ambient,  whereas  at  770K,  p  continues  to  increase  for  an 
interval  following  irradiation,  and  after  several  hundred  microseconds 
levels  off  at  a  value  some  9%  above  the  original  value.  The  latter  effect  is 
ascribed  to  an  increase  in  structural  defect  density  during  the  high  temperature 
anneal. 


1 
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Belozerov,  S.  A.,  G.  M.  Zverev,  V,  S.  Naumov, 
and  V,  S,  Pashkov,  Destruction  of  transparent 
dielectrics  by  radiatljn  from  a  mode-locked  laser. 
ZhETF,  V.  62,  no.  i,~1972,  294-299. 


A  comparative  test  was  made  on  damage  thresholds  in  various 
dielectrics  from  pulsed  laser  radiation.  The  main  purpose  was  to  illustrate 
the  difference  in  effect  of  single  pulse  (10  ns)  and  a  train  of  short  pulses 
(30  X  4  ns)  from  a  mode-locked  Nd  laser  on  dielectric  breakdown.  The 
test  materials  included  type  K-8  glass,  fused  and  crystal  quartz,  leucosapphi 
and  ruby  with  and  without  color  centers.  The  test  conditions  and  findings 
generally  duplicate  those  of  Orlov  et  al  (Effects  of  High  Power  Lasers. 

Dec.  1971,  p.  49),  who  also  advised  on  the  present  experiment.  A  typical 
ilamentary  breakdown  in  the  short-pulse  regime  indicated  self- focusing, 
evidently  not  thermal,  which  was  not  apparent  in  the  monopulse  regime. 

e  tests  verify  that  while  energy  densities  for  threshold  are  of  the  same 
order  of  magnitude  for  all  dielectrics  tested,  the  power  density  threshold 
for  the  short- pulse  mode  is  typically  several  orders  higher  than  for  the 
10  ns  mohopulse,  and  was  found  to  go  as  high  as  lO^^  This  suggests 

a  thermal  relief  mechanism  operating  between  pulses  in  the  pulse  train 
case.  Photos  are  also  included  comparing  filament  appearance  of  the  ruby 
with  color  centers  to  that  in  the  remaining  specimens. 


re. 


Zverev,  G.  M. ,  Ye.  A.  Levchuk,  V.  A.  Pashkov, 
and  Yu.  D.  Poryedin.  Optical  destruction  of  the 
surface  of  lithium  niobate.  ZhETF  v  A?  i 
1972,  307-312. - 


Anomalous  breakdown  thresholds  of  laser- irradiated  LiNbO, 
are  exarnmed  and  the  types  of  destruction  mechanisms  taking  place  are 

lowf^h  =°"t^st  to  most  dielectrics,  LiNbCs  has  both  a  markedly 

ower  breakdown  threshold  at  room  temperature,  as  well  as  a  distinctly 

\  ^".breakdown  level  with  increase  in  ambient  temperatLe. 

User  r/nr  ^^diation  tests  with  a  focused  Q-switched  Nd  glass 

laser  at  1.06^on  polished  LiNbOs  specimens,  in  which  20  ns,  0.  1  j  pulses 

increa^^e^  ^  temperature  (121)  Mw/cm2) 

at  the  with  ambient  temperature  as  seen  in  Fig.  1,  exhibiting  step  jumps 
urie  points.  In  contrast,  the  threshold  characteristic  of  LiTa03 
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Fig,  1,  Relative  threshold  intensity  vs, 
temperature,  LiNb03, 

and  BaTi03  were  found  to  have  no  temperature  dependence  until  the 
Curie  point  was  reached,  while  for  ruby  and  glass,,  no  temperature 
dependence  of  threshold  could  be  found  between  20--700OC,  At  lower 
temperatures  a  cumulative  effect  of  laser  pulses  on  threshold  is  also 
noted  in  LiNbOa  which  is  shown  in  Fig.  2.  An  analysis  of  these  findings 


Fig.  2.  Relative  threshold  vs.  total  laser 
pulses,  N, 

indicates  that  for  temperatures  below  330OC,  surface  damage  is  a  function 
of  light  absorption  in  small  trapping  centers,  whereas  at  higher  temperatures 
the  predominant  mechanism  is  absorption  by  free  carriers  in  the  surface 
layer. 
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Uglov,  A.  A.  Work  presented  in  the  seminar 
on  the  physics  and  chemistry  of  materials 
processing  by  concentrated  energy  beamw- 
FlKhOM,  no.  5,  1971,  158-159. 

on  the  selected  papers  In  the  27th  Seminar 

le  subject,  held  m  Moscow  during  February  1971  Tn*  am; 

wa.  chaired  by  Academician  N.  N.  Bybaii'.  and  at.^ded  b;  ovlr  11 
leading  researcher,  in  the  field  from  Soviet  institutes.  The  bulk  of  the 

CO  7a  ‘’‘f 

nJrf  1  K  ^  mentioned.  To  a  lesser  degree 

particle-beam  interactions  are  also  discussed.  ^ 

oxide  filmf  ^7  y*  discussed  the  growth  kinetics  of  thin 

film  rn  f  1  ®  heating,  and  methods  for  optimizing 

ilm  control.  Experiments  with  a  c-w  CO^  laser  have  yielded  a  40-50^A 

caltla.  ^  -  chrome  at  1  millisec  exposure,  which  is  consistent  with 
calculated  values.  A  number  of  test  results  by  the  same  authors  on  polv- 
mer  targets  was  also  discussed.  ”  ' 

of  dielectrfcs  bv  presented  their  findings  on  vaporization 

ontiral  ^lll?  H  q'^artz  and  other 

have  devl!oLd"f®'  ^  ^  ^^^Posure.  From  their  results  the  authors 

have  de^eloped  an  approximate  mathematical  model  based  on  the  law  of 

mass  conserva  ion  Commenting  on  this,  Yu.  N.  Lokhov  noted  thrsimilarity 
in  characteristics  between  CO2  destruction  of  dielectrics  and  the  self-^on 
sistent  damage  mode  of  metals"  under  giant  pulse  exposure. 

milliner  fpnt  al  discussed  aspects  of  laser 

milling  techniques  in  microelectronics;  the  authors  have  developed  a 

technique  for  producing  laser-cut  thin  film  resistors. 

An  interesting  report  by  M.  S.  Baranov  et  al  described  computer 
and  experimental  studies  on  f  »  effect  nf  i=.oo-  1  ,  .  computer 

tests  Reoiiifo  s.sjf-w,  i.  .  ”  ®7^  laser  pulse  envelope  in  beam-target 

waveforms  sh  H  ^ec  angular,  triangular  and  approximately  sinusoidal 
Tnd  r!  showed  maximum  metal  vaporization  for  the  triangular  waveform 
and  for  the  half-tvave  sinusoid.  A  criticism  offered®  here  by 

A.  A.  Uglov,  which  applied  generally  to  the  laser  beam-target  tests  was 
a  attention  tends  to  focus  only  on  the  strictly  thermal  effects  occurring. 
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wherc&s  chemico«thermal  mechanisms  also  play  a  part  and  warrant 
more  detailed  examination. 

Techniques  for  electron  beam  processing  of  materials  are 
mentioned  by  A.  N.  Kabanov,  D.  B.  Zvorykin  et  al,  and  Vu.  D. 
Belotovskiy  et  al.  These  concerned  nonthermal  beam  processing, 
electron  lithography,  and  beam  deposition  techniques  for  metal  films. 

In  summary,  “Dr.  Rykalin  observed  that  advances  in  laser 
tunability  are  needed  to  expand  laser  processing  techniques,  and  that 
more  detailed  study  of  chemico-thermal  processes  is  needed  for  both 
laser  and  particle  beam  techniques. 


Anisimov,  S.  I.  ,  and  V.  I.  Fisher.  Ionization 
relaxation  and  light  absorption  behind  a  strong 
shock  wave  in  hydrogen.  ZhTF,  no.  12,  1971, 

2571-2576. 

A  simplified  analysis  is  presented  which  describes  the  effect 
of  ionization  kinetics  behind  a  shock  wave  on  e-m  ladiation  absorption 
in  this  region.  The  model  assumes  a  plane  stationary  shock  wave,  and 
limits  the  consideration  to  hydrogen,  since  this  avoids  the  complications 

of  multiple  ionization.  The  ionization  relaxation  zone  in  the  shock  wave 
wake  is  divided  for  analytical  purposes  into  two  regions  with  differing 
mechanisms  of  free  electron  formation:  (1)  a  region  of  "seed  electron" 
formation  from  atom-atom  collisions;  and  (2)  an  electron  avalanche 
region,  with  electrons  freed  predominantly  by  electron-atom  collisions. 

Most  of  the  optical  absorption,  as  well  as  the  significant  change  in  gas 
state,  occurs  in  region  (2),  which  is  therefore  the  region  mainly  treated 
in  the  paper.  Equations  for  incident  flux  density  and  degree  of  ionization 
are  numerically  integrated,  together  with  equations  for  one-dimensiona] 
stationary  flow  of  the  gas.  The  calculated  results  of  density  and  tempera¬ 
ture  profile,  together  with  ionization  characteristics,  are  shown  to  have 
a  definite  correlation  with  incident  flux  density.  The  authors  emphasize 
that  the  characteristics  of  hydrodynamic  variables  in  such  a  shock  wave, 
i.  e.  one  absorbing  an  intense  optical  flux,  will  differ  from  the  corresponding 
parameters  5n  the  usual  detonation  wave,  owing  to  the  ionization  relaxation 
region  which  will  exist  in  the  former  case.  Some  graphical  solutions  of 
the  results  are  included. 
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Norlnikiy,  L.  V.  Initiation  of  a  controlled 
^«akdown  in  gaa  bTThird-harmonirT^^on 
.a^eodymium  laser.  IN:  Kvantovaya  elektronlka.  • 

1971,"  l08-*ior''*''  *-adio,  no.  5, 

by  Ak.anov"er:r(Z^"L\" 

breakdown  in  cas  vas  triffa».'rJ  k’  i^'  *  directional 

:rcfin‘"at T'r"  'tU  p:::e^ranthi:^:si;iTca"^^^^^^^^^ 

har  I  air,  obtaining  a  controlled  breakdown  from  the  third 

to  repeat  the  effect  at  the  fundame  t  i  rT*  optical  breakdown.  Attempts 

since^reakdow  thr^  —  unsuccessful. 

..  .  niiuiu  was  Iirst  exceeded  in  both  cases.  The  results 

fundamental  (1.  17  ev)  and  second  harmonic 
(2.34  eO  to  generate  .he  controlled  breakdown  condition  owing  to  their 
insufficient  photoionization  levels.  Norinskiy  emphasizes  thft  his  results 

Tnote  mod"  should  be  repeated  with  a 

slngle-mode  regime  for  better  clarification  of  this  phenomenon. 


Kaytmazov,  S.  D. ,  A,  A.  Medvedev,  and  A,  M. 

Prokhorov.  Effect  oPa  400  koe  magnetic  field 
on  the  plasma  of  a  laser  spark.  ZhETF  P  v.  14 
1971,  314-316. 

:  *  ji  j  experiment  is  briefly  described  in  which  the  controlling  effect 

is  studied  of  an  external  magnetic  field  on  the  geometry  of  a  laser  spark 
plasma.  Two  conditions  must  evidently  be  met  for  field  control  of  spark 

geometry,  namely  (1)  field  pressure  must  exceed  gas  kinetic  pressure 
in  the  plasma,  and  (2)  the  skin  layer  should  not  exceed  spark  radius  r 
This  means  that  the  external  field  must  be  sufficiently  great  that  on 
lowering  of  plasma  pressure  to  the  magnetic  pressure  level,  plasma 
temperature  still  remains  high  enough  to  preclude  diffusion  in  the  external 
•fielct.  The  corresponding  threshold  for  field  control  in  the  present  case 
was  calculated  to  be  on  the  order  of  300  koe.  Tests  to  corroborate  this 
Jevels  up  to  500  koe,  using  a  transformer-fed  one-turn  coil 
of  0.  8  cm  dia.  instead  of  the  usual  capacitor  bank.  A  lOO/tsec  field  pulse 
was  thus  generated,  which  simplified  the  requirement  of  exact  synchro¬ 
nization  of  laser  spark  and  field  pulse.  Tests  were  run  in  ambient  air 
using  a  neodymium  glass  laser  at  2-3  j  in  both  giant  pulse  and  spike 
regimes  to  produce  breakdown.  The  comparative  effect  of  the  field  is 
seen  in  Pig.  1,  where  the  spark  is  confined  to  a  cylindrical  form  with 
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a 


Fig.  1.  Field  effect  on  laser 
spark,  a.c  -  no  applied  field; 
b  -  field  applied 


a  smooth  boundary.  In  both  laser  regimes  the  field  increased  spark 
axial  length  by  about  1.  5  times;  it  follows  that  this  formation  should 
retard  plasma  cooling.  Nominal  spark  parameters  of  r  =  0. 1  cm  and  time 
constant  r  =  3  x  lO""^  sec  led  to  the  conclusion  that  the  plasma  temperature 
attained  was  at  least  6  x  10^  deg.  K. 


Afanas'yev,  Yu.  V.,  and  V.  B.  Rozanov. 

Spectrum  of  multiply-charged  ions  in  a  laser 
plasma.  ZhETF,  v.  6Z,  no.  1,  I97Z,  Z47-Z5Z. 

A  physical  model  is  proposed  for  the  energy  spectral  form  of 
multiply  ionized  atoms  in  a  laser  plasma.  Based  on  certain  assumptions 
regarding  plasma  diffusion,  the  model  permits  the  development  of  an 
analytical  expression  which  describes  the  desired  ion  energy  distribution. 
The  analysis  demonstrates  that  a  principal  factor  governing  the  energy 
spectrum  for  a  range  of  Z-charged  ions  is  the  recombination  process 
during  the  diffusion  period  following  termination  of  the  laser  pulse.  It 
is  furthermore,  shown  that  the  set  of  Z  present  at  the  end  of  the  pulse 
does  not  necessarily  contain  all  Z  values  present  following  this  time.  For 
simplicity  a  spherically  symmetrical  plasma  flare  is  assumed  to  exist 
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law  Th.  I  ‘"‘o  “  according  lo  a  solf-.lmllar 

danlljN  •P®'*'"”  *•  a.  a  function  of  inWal  flare 

t  -  I*nn  °  a  '’*>“*•  =  1020/cm3,  temperature 

ttae  ,  ~Vx  lo-’o  *  characteristic  photorecombination 

time  tpj.  3  X  10  ^0  sec,  which  is  less  than  plasma  diffusion  time,  hence 

worrofT  W  It  'f  TT"  appreciable.  Analogoua  expert^nta! 

1971.  1306).  but  lack  of  complete  data  from  the  latter  preclude  a  useful 
comparison  of  theory  with  experiment.  A  more  recent  similar  tork  o^ 

if  an  u^nu’”  (P'aaaaa  Phycics,  v.  13,  1971,  19)  in  which  decay 

of  an  LiH  plasma  was  calculated  from  ionisation  and  recombination  processes. 

Yevtushenko,  T.  P. ,  v.  Kh.  Mkrtchyan,  and 
.  .  Ostrovskaya,  Spectroscopic  stiidie.  a 

^ser  spark.  IV.  Absorplion  spectrum  of  a  snei’k 
ItLhytirogen.  ZhTF,  no.  12,  1971,  2581-2589. 

in - - - P*  7,  on  spark  spectra  in  air  hL  a  ' — ^ 


Fig.  1.  Schematic  for  spectral  study  of  a 
laser  spark  in  hydrogen.  Aj.Aj  .  monitor  tubes; 
C  .  hydrogen  vessel  at  6  atm;  F,.  Fj  .  sparks; 
Si,  52  -  spectrographs. 


o 
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and  develapcd  giant  pulses  of  0.  5  j  using  a  driven  prlem  Q- switch 
rotating  at  27, 000  rpm.  A  large  amount  of  graphical  data  on  spark 
paramoterii,  deduced  from  the  spectral  response,  is  presented  and 
analyzed  including  the  time  characteristics  of  absorption,  N.,  and 
plasma  temperature.  An  example  is  given  in  Fig.  2,  showing  temperature, 
pressure  and  density  variation  in  the  center  of  the  spark.  In  conclusion 


Fig.  2.  T,  p,  p  (t)  in  the  center  of  a  laser 
spark  in  hydrogen.  1  -  pressure;  2  -  density; 

3  -  temper’  iure. 

the  authors  discuss  the  accuracy  of  the  spectral  method  for  deducing 
the,  cited  spark  parameters. 
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XIII.  MECHANO-ELECTROMAGNETIC  EFFECT 


Tyarikova,  L.  A. ,  B.  G.  Averbukh,  N.  I.  Mor.kvitin, 
and  N.  A,  Krotova.  Study  of  the  parameters  of  r-f 
radiation  from  breakdown  of  a  polymer  bond  to  a 
solid,  dan  SSSR,  v.  201,  no.  4,  1971,  833-836. 

The  relationship  of  adhesion  mechanical  properties  to  r-f  radiation 
parameters  in  film  +  steel  and  film  +  glass  bonds  was  studied  under  atmospheric 
conditions.  The  films  were:  KLT  polymer  (a  mixture  of  polyethylene  terephthalate 
and  natural  rubbers),  adhesive  plaster,  cellulose  esters,  knd  PVC,  An  experimental 
apparatus  was  used  to:  (1)  strip  the  polymer  film  at  a  variable  rate,  (2)  measure 
the  energy  of  adhesion  and  duration  of  film  stripping  as  well  as  emitted  r-f  pulse 
length;  (3)  photo  record  radio  signals;  (4)  determine  emission  frequency;  and 
(5)  study  the  shape  and  estimate  the  length  of  ^he  r-f  pulses.  The  radiation  detector 
(Fig.  1)  permitted  observation  of  pulses  with  a  rise  time  of  at  least  0. 1  sec. 


Fig.  1.  Diagram  of  radio  wave  recorder; 

1  -  antenna,  2  -  input  circuits,  3  -  mixer, 

4  -  heterodyne,  5  -  I-F,  6  -  detector, 

7  -  oscilloscope. 

Sensitivity  was  30/UV ,  transmission  band  width  was  Af  =  10  kHz,  and  the 
intermediate  frequency  was  f  =  460  kHz.  A  1  m  telescopic  antenna  was  placed 
at  a  maximum  distance  of  1,5  m  from  the  film  separation  opening. 

Pulsed  radiation  of  a  minimum  1  MHz  frequency  was  detected  during 
stripping  of  all  of  the  films,  with  the  exception  of  cellulose  esters.  A  pulse 
was^detected  however,  during  the  instantaneous  separation  of  a  cellulose  aceto- 
butyrate  film  from  glass.  The  rate  of  stripping  consequently  is  a  definite  factor 
in  the  generation  of  r-f.  The  intensity  I  of  radiation  increases  with  an  increase 
of  dxscharge  current  from  the  electrical  double  layer.  This  was  confirmed  by  an 
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Fig.  2.  Energy  A  (1)  and  radiation 
intensity  I  (2)  versus  peeling  rate  v  of 
a  KLT  rubber  film  from  a  metal. 

experimentally  observed  increase  in  pulse  amplitude  with  increases  in  stripping 
rate  v.  At  a  given  v,  I  is  higher  when  adhesion  is  stronger  (Fig.  2).  The 
magnitude  of  I  in  film  separation  from  metal  is  lower  than  in  dielectric  separation. 
An  oscilloscope  trace  of  10  MHz  radiation  during  stripping  of  a  KLT  film  from 
a  dielectric,  e.  g.  glacs,  at  v  =  I  m/sec  shows  high-amplitude  periodic  pulse 
packets  in  addition  to  smaller  amplitude  pulses.  The  period  T  between  high 
amplitude  pulses  is  constant,  i.  e.  dynamic  resistance  is  constant  for  a  given 
gas  and  gas  pressure.  Electromagnetic  phenomena  observed  in  adhesion 
breakdown  generally  follow  the  same  pattern  as  in  a  gas  discharge,  but  also 
display  certain  unique  characteristics.  The  presence  of  the  interval  T  is  related 
to  accumulation  of  a  static  charge  along  a  considerable  portion  of  the  stripped 
film;  surface  electrical  conductivity  must  thus  play  an  important  part  in 
adhesion  breakdown.  The  time  dependence  of  radiation  amplitude  remains  uniform 
within  a  1-10  MHz  range  of  frequencies.  Radiation  intensity  oscillations  in  sectors 
of  small  and  high-amplitude  pulses  re’/^eal  the  discrete  nature  of  discharges 
which  accur  in  microregions  of  a  stripped  film,  where  field  intensity  is  critical. 

In  agreement  with  earlier  findings,  the  length  of  an  elementary 
stripped  region  of  rubber  film  on  glass  was  determined  to  be  30^^..  Generation  of 
r-f  oscillations  is  considered  to  be  the  most  convincing  proof  of  gas  discharge 
development  in  a  propagating  crack.  Visible  radiation  and  acoustic  vibrations 
were  recorded  along  with  r-f  in  all  cases.  At  stripping  rates  of  0.  1  -  1  cm/sec, 
the  visible  and  r-f  pulses  were  synchronized  as  determined  from  an  oscilloscope 
trace  of  rubber  film  stripping  at  v  =  0.  1  cm/sec.  Oscillatory  emission  of 
"mechanoelectrons"  was  al-j  observed  during  adhesion  breakdown  in  vacuum. 

This  confirms  the  oscillatory  nature  of  the  breakdown  and  associated  phenomena. 
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Abramova,  K,  B,  ,  V.  P,  Valitskiy,  N,  A,  Zlatin, 

B,  P,  Peregud,  and  I.  Ya,  Pukhonto,  Radiation 
generated  by  rapid  deformation  and  rupture  of  metals. 

DAN  SSSR,  V.  201,  no.  6,  1971,  1322-1325. 

Experiments  are  described  which  were  designed  to  test  hypotheses 
on  the  magnetohydrodynamic  origin  of  radiation  emitted  by  rapid  mechanical 
deforn.ation  and  rupture  of  metals.  The  particular  case  considered  was  for  an 
impact  velocity  Vwl,000  m/sec  and  in  the  absence  of  an  electric  field.  Radiation  ir 
the  3,000  -  8,200  A  range  was  not  reliably  detected  in  preliminary  measurements 
of  rupture  of  Cu,  Al,  Mo,  and  Bi  samples  using  an  impact  testing  machine  at 
V  =  7  m/sec.  Other  tests  at  V  up  to  1,200  m/sec  were  done  by  the  technique  of 
"contactless"  rupture  of  a  san^ple  (1)  separating  two  cylindrical  chambers  (2) 
and  (3)  as  shown  in  Fig.  la.  Pressure  was  maintained  at  1  torr  in  chamber  (2) 
and  at  10“^  -  10”2  torr  in  chamber  (3).  Deformation  or  rupture  of  the  rear  side 
surface  (1)  which  closes  chamber  (3)  was  initiated  by  a  15  mm  dia  copper 
impactor  (4).  Radiation  in  the  3,  000  -  8,200  A  range  was  recorded  by  a  photo¬ 
multiplier  (6)  through  window  (7).  A  copper  deformation  rate  of  about  100  m/sec 
was  deterniined  from  x-ray  shadow  photographs  taken  at  regular  intervals. 
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Fig.  1.  Experimental  configurations. 

Rupture  of  the  rear  surface  was  accomplished  in  112/6sec.  At  a  1200  m/sec 
impact  velocity,  oscilloscope  signal  traces  from  the  photomultiplier  show 
2-2.5  msec  pulses  with  several  characteristic  peaks.  Pulse  intensity  from  a 
ruptured  copper  surface  was  3-8  times  that  from  a  deformed  surface  without 
any  visible  rupture.  The  radiation  pulses  were  recorded  with  27  and  33  mm 
thick  samples.  In  one  control  experiment,  a  50  mm  thick  sample  was  not  deformed 
and  no  radiation  was  detected.  A  special  control  experiment  was  set  up  involving 
exclusion  of  triboluminescence  of  the  oxide  film.  In  this  experiment,  a  supposedly 
oxide-free  copper  surface  formed  on  the  rear  face  of  a  27  mm  thick  sample  by 
a  previous  impact  was  reimpacted  2-3  min  after  the  first  impact.  A  light  pulse 
of  the  same  intensity,  but  tixiorter  than  the  first  one,  was  recorded  in  this 
experiment.  In  both  experiments,  radiation  from  collision  or  friction  of  the 
sample  with  the  small  chamber  walls  did  not  interfere  with  the  recordings, 
because  the  colliding  and  friction  parts  were  shielded  from  the  photomultiplier. 
Residual  gas  Interference  in  the  shock  wave  front  was  also  excluded,  since  a  shock 
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wave  could  not  develop  at  subsonic  deformation  velocities.  No  trace  of  crack, 
edge  melting  was  detected  on  the  photographs,  indicating  that  thermal  radiation 
was  recorded.  Analogous  results  were  obtained  in  experiments  with  duralumin. 
Final  proof  of  radiation  emission  from  the  part  of  the  sample  in  the  process  of 
eformation  and  rupture  was  obtained  by  photographing  the  self- radiation  of  a 
luminescent  crack  in  copper  using  an  image  converter  (Fig.  lb).  The  radiation 
contour  in  the  photograph  coincides  closely  with  the  crack  contour. 

u  j  r  condtuded  that  the  experimentally  observed  radiation  is  emitted 

by  deformation  and  rupture  of  the  metal.  This  radiation  exhibits  the  principal 
characteristics  of  luminescence,  i.  e.  it  exceeds  thermal  radiation  at  a  given 
temperature,  and  its  lifetime  exceeds  significantly  the  period  of  a  light  wave. 
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XIV.  MISCELLANF.fVtIR 


Poltavtsev,  Yu,  G, ,  V,  P,  Zakharov -and 
V,  N.  Chugayev.  Structural  studies  of  graphitization 
of  thin  carbon  films  induced  by  powerful  light 
pulses,  Krlstallogira^lya^  v.  16,  no,  2,  1971, 

415-419. 

An  electronographic  analysis  of  graphitization  kinetics  in  thin 
amorphous  carbon  films  exposed  to  powerful  light  pulses  was  performed.  The 

optical  source  in  this  case  was  an  IFP  type  xenon  flashlamp,  operating  at 
voltages  up  to  1.7  kv.  Carbon  films  were  obtained  by  thermal  diffusion  of 
carbon  rods  onto  cold  glass  substrates;  samples  were  then  placed  on  a  copper 
grid  and  exposed  to  various  energy  light  pulses,  and  in  this  manner  various 
stages  of  graphitization  were  obtained.  Kinetics  of  graphitization  was  studied 
from  the  nearest  order  parameters,  which  were  determined  for  each  graphi¬ 
tization  stage  from  the  radial  atomic  distribution  curves.  Radial  atomic 
distribution  curves  were  calculated  from  the  electronograms  obtained  by 
means  of  an  electron  microscope  and  processed  according  to  techniques 
described  by  I.  D.  Nabitovich  et  al.  (Kristallografiya,  no.  4,  1967,  584.) 
Experimental  data  and  theoretical  calculations  indicate  that  there  is  a  definite 
similarity  in  the  structure  of  amorphous  and  graphitized  carbon  films. 
Amorphous  carbon  has  a  partially  ordered  hexagonal  lattice  structure,  similar 
to  that  of  the  graphite  cryst^line  lattice.  However,  there  are  different 
atomic  layer  arrangements  in  amorphous  carbon  and  graphite.  Based  on  the 
analysis  of  the  structure  of  original  amorphous  carbon  films  and  radial  atomic 
distribution  curves,  the  authors  suggest  the  following  graphitization  kinetics. 
During  the  first  graphitization  stages  a  certain  number  of  carbon  hexagons 
with  a  single  C-C  bond  are  destroyed.  Predominantly  double  and  triple* 
carbon  bonds  are  formed  among  the  released  carbon  atoms.  During  the 
subsequent  graphitization  stages  the  hexAgonal  and  structural  ordering  takes 
place.  The  atomic  layers  are  formed  in  such  a  manner  that  every  atom  has 
one  neighbor  in  the  adjacent  layers.  (In  amorphous  stage  every  atom  has  on 
the  average  four  neighbors  in  adjacent  layers.  )  In  the  final  graphitization 
stages  the  hexagons  get  somewhat  distorted;  however,  there  is  no  change  in 
relative  positions  of  the  atoms.  It  is  suggested  that  the  cause  for  the  hexagonal 
distortion  is  the  formation  of  double  carbon  bonds  between  certain  adjacent 
carbon  atoms,  which  thus  coexist  with  single  carbon  bonds  in  the  final 
graphitization  stage. 


77 


I? 


Zaytsev,  A,  V,  A  new  theory  of  ball 
lightning.  ZhTF,  no.  1,  1972,  213-216. 

The  author  expcmnJb  a  theory  of  ball  lightning  formation  based 
on  published  research  data  on  the  formation  of -extremely  fine  electroconductive 
structures  in  an  electric  field.  Such  structures  are  composed  of  fine  metallic 
CT  partly  trtd  bocf  paiW-ele^;.  Analugoira  struotaxea  oombined  with 

a  corona  space  charge  may  form  a  streak  lightning  embryo.  Discharge  current 
from  the  structure  in  a  static  earth  -  cloud  electric  field  flows  through  the 
structure  resistance  and  heats  it.  The  heated  structure  then  emits  visible 
light  which  might  be  perceived  as  ball  lightning.  The  embryonic  structure 
is  visualized  in  the  form  of  a  nucleus  having  a  spatial  network  structure  with 
spatially  interconnected  chain  structures  grown  out  of  it  (see  Fig.  1).  The 


Fig.  1.  Spatial  distribution  of  ball  lightning  structures: 

1  -  spatial  netv/ork  structure  and  space  charge;  2  -  upper 
structures  and  space  charge;  3  -  lower  structures  and 
space  charge;  4  -  separate  ascending  lower  structure 
chains;  5  -  lower  discharge  boundary  of  lower  structures; 

6  -  upper  discharge  boundary  of  upper  structures  in  the 
pre-leader  discharge  phase;  7  -  leader  channel;  8  -  charges 
ascending  to  the  cloud. 
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mechanism  of  formation  of  a  hill  lightning  is  described  on  the  basis  of 
this  model  as  follows.  Discharge  currents  from  separate  chainrof  the 
upper  and  lower  structures  accumulate  on  the  nucleus.  Given  a  sufficient 
number  of  chains,  of  the  appropriate  length,  the  total  current  across  a 

L"  sTon"  A  ^th°  '"T  structure  to  the  point  of  visible  light 

observed*  A  f  ft  ^  luminescence  of  the  nucleus  may  be 

observed.  A  further  increase  in  current  across  the  upper  and  loLr 

structures  results  in  the  accumulation  of  a  significant  space  charge  within 
th^e  embryonic  structures,  which  periodically  forms  a  sLrt  curre^nt  leader 
channel  of  streak  lightning.  The  leader  current  flows  across  tL^mbtyonL 

flashes  and  explosions  of  the  luminous  nucleus.  A  long  leader,  rcsultine 
from  a  sudden  decrease  in  structural  resistance  at  the^nstant  of "art^l" 
of  the"f’  perceived  as  a  strong  explosion  along  with  the  disappearance 
of  the  luminous  nucleus.  The  latter  appears  to  be  the  ball  lightning  energy 

TheT’d^  earth-cloud  electric  field. 

tht  struclurrs'^^Ch  ^  distance  from 

fall  and  f  embryonic  structures  sometimes  alternately 

ilj  f  f/  u  ^  "jumping  ball-  pattern  without  breaking 

7  embryo  moves  only  along  the  earth's 

surface,  subject  to  the  mutually  balancing  forces  of  weight  and  charge 
ttraction  by  the  cloud.  In  this  case  the  charge  of  the  lower  structuL 
breaking  contact  with  the  earth  is  neutralized  by  the  discharge  current- 

pml  ’  Predominance  of  the  charge  attraction  over  the  weight  of  the 

embryomc  structure  manifests  itself  if  the  structure  is  heated  sufficiently 

u current  to  rise  in  ascending  air  streams.  In  that 
case,  the  embryo  may  be  lifted  far  from  the  earth. 

f  Optimum  conditions  for  ball  lightning  formation  are  a 
ombination  of  high  particle  concentration,  a  strong  electric  field  and 
igh  temperature.  Such  a  combination  may  occur  when  a  layer  of  ash 
or  soot  particles  is  polarized  by  a  discharge  (10^  .  iQ?  v)  of  streak 
lightning.  The  polarized  particles  form  a  nucleus  which  is  sintered  into 

discharges.  The  existence  of  a 

y  particles  IS  however  not  a  necessary  precondition.  The  particles 

may  form  in  an  induced  discharge  channel  by  partial  combustion  or  vapor¬ 
ization  ot  carbon-containing  material  or  metals.  Growth  of  the  upper 
structures  in  this  case  is  explained  by  accumulation  after  a  discharge  of 
the  residual  particles  on  the  polarized  upper  part  of  a  spatial  structLe. 
ong  e  ectric  field  created  by  a  current  breakdown  on  the  attached 

growth  of  the  structures.  The  burned-out  spatial  structures  are  recon  - 
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structed  from  ground  material  particles  by  current  breaks  generated 
between  the  spatial  structure  and  earth  at  the  instant  of  separation  from 
the  earth.  The  particles  remaining  below  the  spatial  structure  after  a 
discharge  form  the  lower  structures  and  close  the  discharge  circuit 
across  the  nucleus  to  the  earth.  Ball  lightning  lifetime  is  limited  only 
by  those  conditions  favorable  to  structure  formation,  as  long  as  the 
nucleus  structure  is  not  burned  out  completely. 

Examples  are  given  of  other  ball  lightning  characteristics 
which  might  be  forecast  on  the  basis  of  the  cited  mechanism. 


Kotenko,  V.  G.  ,  and  Yu.  A.  Litvinenko.  Generation 
of  strong  pulsed  magnetic  fields  in  solenoids  of 
bimetallic  wires.  ZhTF,  no.  1,  1972,  183-186. 

The  possibility  of  producing  a  strong  pulsed  magnetic  field  in 
solenoids  has  been  investigated  by  using  reinforced  bimetallic  wire  as  the 
winding  material.  One  bimetal  component  forms  a  highly  electrically  conductive 
part,  with  the  other  being  highly  rigid,  thus  forming  a  reinforced  winding. 
Experience  in  using  such  a  winding  is  outlined.  The  main  design  problem  is  to 
transfer  electrodynamic  stress  produced  in  the  conducting  element  to  the  more 
rigid  part.  The  degree  of  loading  of  the  rigid  part  depends  on  the  construction 
and  the  bonding  characteristics  between  the  two  component  metals.  Theoretical 
relationships  are  derived  for  determininjj  the  amplitude  of  magnetic  field 
corresponding  to  the  stability  limit,  and  for  the  tolerable  pulse  duration  of  the 
magnetic  field.  Computed  relationships  are  reduced  to  the  calculations  of  a 
nonreinforced  ordinary  solenoid,  i.  e.  ,  using  only  the  current-carrying  part  of 
the  bimetal,  sb  that  the  basic  principle  can  be  applied  to  any  type  of  reinforcing 
element.  Fig.  1  shows  the  form  of  characteristic  curves  for  the  relationships 
obtained. 


Fig.  1.  Design  limits  of  bimetal  winding. 

Hi  2  -  magnetic  field  corresponding  to  the 
stability  limit  of  the  reinforced  solenoid; 

Hi  -  magnetic  field  corresponding  to  the  yield 
point  of  ordinary  nonreinforced  solenoid; 

X  =  S2/S1*  where  83  -  cross-sectional  area  of 
rigid  metal  and  Sj  -  cross-sectional  area  of 
conducting  metal;  t  1  2  "  duration  of 

magnetic  field  Hj  2  ^1  "  pulse  duration  of 

magnetic  field  Hj 

^1,2=  ^1 
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Experiments  were  conducted  with  a  NbZrTi  (main)  -  'copper 
(coating)  wire,  having  stability  limits  <^21im'®  ^00  kg/mm^  and  <ri«20  kg/mm^, 
respectively.  At  x  =  1,  the  magnetic  field  can  be  increased  by  2.  3  times,  while 
pulse  duration  decreases  by  20  times.  Results  of  experiments  using  liquid 
nitrogen  and  liquid  helium  along  with  the  calculated  results  are  tabulated  (Table  1). 
The  experimental  results  are  in  good  agreement  with  the  calculated  ones,  proving 
the  applicability  of  bimetallic  wires  for  winding  materials  of  a  reinforced  solenoid. 


Cooling 

Calcu. 

lation 

'i.2 

Expt 

»t,2 

Calcu- 

lation 

Expt 

No.  of 
Pulses 

Remarks 

N 

0.3 

0.25 

460 

410 

-  15 

^Thermal 

. 

.....  - 

_ 

/destruction  of 

He  • 

0.55 

0.25 

460 

500 

~10 

/wires. 

Table  1. 

*  cooling  of  solenoids  lower  than  liquid  hydrogen 
is  not  considered. 

**  voltage  of  capacitor  bank  -  4  kv. 

From  the  experimental  and  calculated  results,  it  may  be  deduced  that  use  of 
aluminum  as  one  and  steel,  tungsten,  etc.  as  the  other  constituent  of  the  bimetallic 
wire  may  make  it  possible  to  obtain  a  relatively  long  pulse  duration  of  magnetic 
field  on  the  order  of  500~600  koe.  A  further  selection  of  optimum  solenoid 
forms  may  increase  the  magnetic  field  by  an  additional  1.2~f.  5  times. 


Alekseyev,  E.  I.  A  description  of  the  interaction 
£^nonc9herent  radiation  with  matter.  ZhPS, 

V.  15,  no.  6,  1971,  1090-1093. 


.  ,  general  probability  expression  is  presented  from  which  the 

interaction  of  radiation  with  matter  can  be  approximately  described  The 
model  eeaumee  moe.  types  of  noncoherent  radiation  sources  of  ma^o 

ranSL"  ■  ‘VP''  “f radiation  as  well; 

rala«^n"°‘r  ‘"’"’ever  excluded.  A  general  expression  defining 

which  c  elements  is  given, 

■.U  treated  as  a  system  of  ordinary  stochastic  differential  equations 

^th  stipulated  or  random  initial  conditions."  A  rigorous  solutron  for  the 
sume  model  could  be  obtained  by  a  form  of  the  Fokker-Planck- 

soMon  ml  P'-»c«cal  considerations  a  satisfactory 

yea  y  using  only  the  mean  value  of  density  matrix  elements. 
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The  author  illustrates  the  approach  by  examples  of  approximate  equations 
describing  differences  of  population  levels  in  a  two-level  system,  under 
the  effect  of  radiation  taken  to  be  a  normal  stationary  process  with  zero 
mean.  The  method  is  tedious,  requiring  computer  solution,  and  is  limited 
to  systems  with  a  few  levels  only,  but  in  theory  applies  to  radiative 
interaction  of  arbitrary  intensity  and  spectral  characteristic  with  matter, 
as  long  as  the  radiation  may  be  considered  as  a  normal  stationary  process. 


Askar'yan,  G.  A.  ,  and  T.  G.  Rakhmanina. 

Scattering, refraction  and  reflection  of  sound 
under  the  action  of  intense  light  on  a  medium. 

ZhETF,  V.  61,  no.  3,  1971,1199-1202. 

An  analytical  discussion  is  given  describing  the  effect  of 
powerful  optical  absorption  in  a  medium  on  acoustic  parameters  of  that 
medium.  Coherent  or  noncoherent  radiation  of  sufficient  intensity  can 
cause  substantial  local  changes  in  acoustical  scattering,  refraction  and 
reflection  owing  to  gas  evolution,  bubble  formation,  etc.  The  case  of  a 
liquid  medium  is  treated  here,  with  bubbles  assumed  to  be  formed  by  beam 
absorption,  as  this  produces  drastic  changes  in  the  cited  acoustical  para¬ 
meters.  An  expression  for  mean  bubble  size  is  derived  as  a  function  of 
breakdown  threshold  and  other  parameters  of  the  liquid  medium.  A  compar¬ 
ison  of  theoretical  scatter  in  bubbles  with  scatter  in  locally  heated  regions, 
at  the  same  per  unit  energy  absorption,  shows  that  the  bubble  scattering 
cross-section  will  exceed  that  of  tlie  heated  regions  by  as  much  as  10^^, 
which  indicates  the  predominant  scattering  effect  of  bubble  formations. 

A  brief  treatment  on  ijound  refraction  and  reflection  from  optically 
disturbed  regions  is  also  given.  The  discussion  in  general  is  concerned 
with  acoustic  velocities  ranging  from  supersonic  to  hypersonic. 


